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Abstract 


In this paper we describe the history of the LHCb experiment over the last three 
decades, and its remarkable successes and achievements. LHCb was conceived 
primarily as a b-physics experiment, dedicated to CP violation studies and measure- 
ments of very rare b decays, however the tremendous potential for c-physics was also 
clear. At first data taking, the versatility of the experiment as a general-purpose 
detector in the forward region also became evident, with measurements achievable 
such as electroweak physics, jets and new particle searches in open states. These 
were facilitated by the excellent capability of the detector to identify muons and to 
reconstruct decay vertices close to the primary pp interaction region. 

By the end of the LHC Run 2 in 2018, before the accelerator paused for its second 
long shut down, LHCb had measured the CKM quark mixing matrix elements and 
CP violation parameters to world-leading precision in the heavy-quark systems. 
The experiment had also measured many rare decays of b and c quark mesons and 
baryons to below their Standard Model expectations, some down to branching ratios 
of order 107°. In addition, world knowledge of b and c spectroscopy had improved 
significantly through discoveries of many new resonances already anticipated in the 
quark model, and also adding new exotic four and five quark states. 

The paper describes the evolution of the LHCb detector, from conception to its 
operation at the present time. The authors’ subjective summary of the experiment’s 
important contributions is then presented, demonstrating the wide domain of 
successful physics measurements that have been achieved over the years. 
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1 Introduction 


LHCb is an experiment at the CERN LHC, dedicated to the study of heavy flavours 
with large statistics. The resulting high precision makes possible the observation of tiny 
deviations from the predictions of the Standard Model (SM) in CP violation and rare 
phenomena, variations which could hint at New Physics (NP) processes. LHCb started 
taking data in 2010. The so-called Run 1 commenced at an initial centre of mass energy of 
Vs = 7 TeV which was then increased to ys = 8 TeV, collecting an integrated luminosity 
of 3.23 fb! until the end of 2012. After a two-year shutdown, LHC operation continued 
from 2015 to 2018 (Run 2), when the experiment took data at \/s = 13 TeV, recording 
an integrated luminosity of ~ 6 fb~'. Throughout the running periods, LHCb collected 
and analysed an unprecedented number of b decays, and also enlarged its scope to include 
charm physics, W and Z measurements, jets and nuclear collisions. 

In this paper the motivation for the LHCb experiment is described, recalling how its 
design was developed and evolved. The experiment’s major achievements in terms of 
physics results are then summarised. Finally we discuss the plans for the future upgrade, 
in the period when the Super-KEKB collider will also operate. 

The layout of this paper is as follows. The Introduction (Sect. 1) describes the 
status of the b physics programme at the time when the LHCb detector was conceived, 
and provides an evolution of its design. In Sect. 2 the basic elements of the detector 
are described, optimised for the requirements of heavy-flavour physics measurements, 
together with a description of the triggers. The measurements which were originally the 
major aims of the experiment, i.e. the CKM matrix measurements and CP violation 
and very rare decays of the b-quark, are reviewed in Sects. 3 and 4, respectively. In 
Sect. 5, a summary of the wide-ranging results in b- and c- spectroscopy is presented. 
A review is given of the many non planned physics areas in Sect. 6, such as results 
on jets, electro-weak (EW) physics, and searches for new particles not associated with 
heavy flavours. In all these domains, LHCb proves to be an extremely versatile detector, 
providing complementary measurements to those of the LHC General Purpose Detectors 
(GPDs). The paper concludes in Sect. 8 with a short description of the upgrade plans, 
which will ensure LHCb operation beyond 2030. 

The paper presents the authors’ subjective summary of LHCb’s many major physics 
results, however the review inevitably omits a substantial number of important mea- 
surements. To this end, additional information can be found in the paper’s exhaustive 
bibliography. 


1.1 b physics at the end of the XXth century 


In 1970 Glashow, Iliopoulos and Maiani postulated the existence of a fourth quark, 
charm, necessary to explain the smallness of K? oscillations in the framework of the 
Standard Model of Weak Interactions (i). The so-called Glashow-Iliopoulos-Maiani (GIM) 
mechanism generalized Cabibbo’s idea of rotated weak currents to the quark model with 
two doublets, introducing a 2 x 2 mixing matrix written in terms of the Cabibbo angle. An 
experimental confirmation of the fourth quark hypothesis came in 1974 with the discovery 
of J/w 2/3}, soon identified as a cé bound state, and followed by the discovery of open 
charm. In 1973 Kobayashi and Maskawa proposed a third heavy-quark doublet in 
order to describe CP violation in the framework of the SM, thus generalizing the Cabibbo 


matrix to the 3 x 3 CKM (Cabibbo, Kobayashi, Maskawa) matrix. The discovery of the 
Y in 1975, followed by the charged and neutral B-mesons in 1983 |5| proved the validity 
of their idea and held the prospect of understanding quantitatively CP violation. 

The level of CP violation in b quark decays was expected to be orders of magnitude 
larger than in the neutral K system, but unfortunately the relevant decay channels had 
only tiny branching fractions, so the lack of intense “b-quark sources” slowed the progress 
of beauty physics: in particular, in 1986 the PDG only listed five decay modes of B° and 
B+. However in 1987 a new impetus came from the discovery at ARGUS of B° — B 
oscillations (6. It was clear that the forthcoming LEP machine, designed for an entirely 
different purpose, and the symmetrical CESR collider, could not yield an exhaustive answer 
to all the questions related to the CKM hypothesis, despite their valuable contributions 
to many facets of b-physics (7/3). When in 1989 P. Oddone proposed an asymmetric 
ete” collider [o] operating at the Y (4S) energy with a luminosity above 1033 cm~*s7!, an 
intense period of accelerator studies ensued. This gave birth to the PEP-II and KEKB B 
Factories, which were approved in 1994, and started operating in 1998, soon reaching and 
passing their design luminosity. 

Around this time, proponents pursued the idea of exploiting hadron beams to attack 
the problem of detecting CP violation in the b sector. The idea behind this was that 
the large hadronic b production cross-section plus the high-intensity hadron beams at 
the already existing and planned proton accelerators would produce a large number of 
bb pairs, sufficient to gather evidence for CP violation at least in the so-called “golden 
channel” B° + J/wK§. To achieve an adequate background rejection, the experimental 
difficulties were formidable because of the small ratio of the b cross-section to the total 
hadronic cross-section at the ys values available in fixed-target and collider experiments. 

In 1985, the fixed-target WA75 hybrid experiment at the CERN SPS observed in 
emulsions the first partially reconstructed bb pair produced by an extracted pion beam 
of 350 GeV/c, confirming that the production cross-section at such low energies was 
very small. To circumvent this problem, simple experiments were proposed for the 
CERN SPS and for the planned UNK machine at Serpukhov, with a brute-force 
approach based on a high-intensity extracted beam and a minimalist detector designed to 
reconstruct the B° — J/WK® decay and to provide flavour tagging (Fig. (1.1). There was 
no time-dependent analysis of the decay since the beam-dump character of the experiments 
made the use of a micro vertex detector impossible. There were exploratory fixed-target 
experiments, at CERN (WA92) and at Fermilab (E653, E672, E771, E789) |14|, which 
tried to observe and measure beauty events, albeit without (or very limited) success. 

In 1989 P. Schlein proposed a dedicated Beauty experiment exploiting the large b 
cross-section expected at the CERN SPS proton-antiproton collider (ys = 630 GeV) (15). 
The authors of the proposal (P238) remarked that the bulk of bb production occurred 
at very small angles with respect to the beams, therefore making a compact experiment 
practical. The heart of the detector was a Silicon Microvertex Detector operating very 
close to the beam (1.5 mm) coupled to a fast readout and track-reconstruction electronics. 
The Microvertex Detector provided the trigger by requiring that accepted events had to be 
inconsistent with a single vertex. P238 was not approved but the CERN R&D Committee, 
established to support new detector developments in view of the LHC, approved in 1991 a 
test of the Microvertex Detector in the SPS Collider. This proved very successful 
and paved the way towards the future COBEX experimental proposal. 

At the time when the ete~ colliders were approved, the HERA-B experiment had been 


Figure 1.1: A sketch of a proposed beam-dump style experiment designed to detect B° — Tpke 
decays. The Ke would traverse the conical slit before decaying [11]. 


conceived at DESY {18]. Approved in 1994, HERA-B exploited the 920 GeV HERA proton 
beam on a fixed target made of metallic wires, placed inside Roman Pots in the vacuum 
pipe, and immersed in the beam halo. HERA-B was approved to take data in 1998, one 
year before PEP-II and KEKB. The sophisticated apparatus consisted of a single-arm 
spectrometer, including a RICH, a large microvertex silicon detector, a high-resolution 
tracker, plus an electromagnetic calorimeter. HERA-B was designed primarily for the 
detection of the B° > J/Y K8} decay and its trigger was based on J/ọY reconstruction at 
the first level. 

At \/s = 40 GeV, the bb cross-section is about 1078 of the total hadronic cross-section, 
hence HERA-B had to achieve a background rejection around 107"! for the B® + J/WK& 
decay. Data taking conditions were similar to those of the current LHCb experiment (a 
40 MHz interaction rate), as were the requirements of radiation resistance. HERA-B 
started data taking in 2000 but it soon emerged that the detector did not have sufficient 
rejection power against background and the track reconstruction was not as efficient as 
expected. The large number of detector stations and their total thickness in terms of 
radiation lengths made secondary interactions an important issue for event reconstruction. 
Eventually HERA-B could not observe b events efficiently, but taught several valuable 
lessons for any future experiment working in a crowded hadron-collision environment. 
These were a need of a robust and efficient tracking and a flexible trigger systems able 
to adapt to harsher environments than may have been expected, as well as the need to 
design the thinnest and lightest detector (in terms of radiation and interaction lengths). 


1.2 Towards the LHC 


Over the same period, the planned LHC and SSC machines, with their large energies, 
promised spectacular increases of the bb cross-section, thus making the task of background 
rejection much simpler. This was particularly true for operation in collider mode, but 
even in fixed-target mode (ys ~ O(130) GeV) the b cross-section was expected to be 


a respectable 1 wb at the LHC and 2 pb at the SSC [19]. The large cross-section and 
corresponding good background rejection would facilitate a hadron B Factory, competitive 
and complementary with e*e~ colliders, which could focus primarily on the measurement 
of CP violation and also allow the study of the spectrum of all b particles. Given the 
intrinsically “democratic” nature of hadronic production, the new hadron machines would 
also give access to large samples of B? and of b baryons, something not possible at ete 
colliders operating at the Y resonances. 

Two schools of thought soon emerged: one pursuing a fixed-target (FT) strategy and 
the other based on a collider mode. The more favourable ratio of the bb to the total 
hadronic cross-section, about two orders of magnitude larger in collider mode, gave this a 
competitive advantage. 

There was, however, a strong reason in favour of the FT concept: a collider B 
experiment could not operate at the design luminosity of the machine (1033734 cm?s~) 
because of the significant number of overlapping interactions (pile-up) with multiple 
vertices. This would have required dedicated low-luminosity running, creating a potential 
conflict with the major experiments and considerably reducing the data-taking time. Later 
it was ascertained that the individual experiment luminosities could be tuned over a broad 
range with an appropriate design of beam optics in the interaction regions, hence this 
would become a moot point, but at the time it was a serious one. 

Moreover, for the advocates of the FT approach, the advantage of the larger cross- 
section in collider mode was partially offset by the higher event multiplicity and by the 
shorter flight path of beauty particles. In addition, while the pr of the b decay products 
would to a good approximation be the same in the two modes, the pr of the other collision 
products would be smaller in FT mode, thus making the trigger simpler. Active silicon 
targets were also possible with an extracted beam, where the track of a charged b-hadron 
would be directly measured. 

Finally, bb production kinematics is forward peaked in the centre-of-mass (CM) system, 
but the Lorentz boost of the centre-of-mass in the FT mode (8y > 60) concentrates the 
event at smaller angles than in collider mode. It was therefore possible, in principle, 
to build a more compact detector, achieving larger angular acceptance at lower cost. 
There was also the possibility of recycling components (in particular dipole magnets) 
from existing detectors. The cost was an important consideration, since a dedicated B 
experiment at the LHC (or SSC) was generally considered to have secondary importance 
with respect to the general-purpose experiments. 

The anticipated demise of the SSC led three groups to study and propose dedicated 
b experiments at the LHC. The three Letters of Intent (Lol) were presented in 1993. 
COBEX [20], an acronym for Collider Beauty Experiment, with P. Schlein as spokesperson, 
was a collider experiment with a backward-forward geometry. The other two proposals, 
Gajet (spokesperson T. Nakada) and LHB (spokesperson G. Carboni) were fixed 
target experiments, the first, as the name suggests, using a gas jet target, the latter 
exploiting an extracted beam. Since no traditional beam extraction was foreseen for the 
LHC, LHB (Large Hadron Beautyfactory) used a parasitic extraction technique, based on 
channeling in a bent silicon crystal placed close to one of the circulating beams (Fig. |1.2). 
A dedicated R&D experiment, RD22, approved by the CERN DRDC to test the 
feasibility of this idea at the SPS, demonstrated that high-efficiency beam extraction 
(larger than 10%) was possible [24]. The three proposed experiments were presented in 
their final form in 1994 at the Beauty’94 Conference [25]. It should also be noted that, 
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Figure 1.2: Proposed LHC beam extraction scheme based on crystal channeling. The channeled 
beam, deflected by 0.7 mrad in a bent silicon crystal, was guided towards the beam transport 
tunnel by several conventional magnets. The cost of the extracted beam was a non-negligible 
fraction of the detector cost. 


by then, CDF had already contributed important b-physics results, offering a glimpse of 
what would later prove to be the extraordinary success of b-physics at hadronic machines. 

Following the submission of the Lols, the LHC Committee (LHCC) considered the 
proposals in June 1994. One of its concerns was the small size of the three collaborations: 
the total number of physicists involved was barely one hundred. In addition the Committee 
remarked that beam extraction by channeling could not be guaranteed at that stage, be- 
cause accelerator experts feared possible interference with normal LHC operation. Finally 
the following recommendation was issued: “The collider model approach has the greater 
potential in view of the very high rate of b production, the much better signal/background 
ratio and the possibility of exploring other physics in the forward direction at 14 TeV”. 
The LHCC encouraged the three collaborations to join together and design a new ex- 
periment, incorporating attributes of each, and operating in collider mode. The LHCC 
noted that its ambitious request was justified by the fact that, at the startup of the LHC, 
the experiment would not be an exploratory one, since CP violation would already have 
been observed “at HERA-B, FNAL or B Factories”. The LHCC also issued guidelines 
requiring a number of issues to be addressed and solved. The three experiments combined 
into the new LHC-B Collaboration (then named), which published its Letter of Intent in 
1995 [26]. 

The new experiment derived several of its characteristics from the parent proposals: 
notably the concept of a silicon vertex detector in retractable pots, the calorimeter design, 
and the high pr first-level trigger. 

In the transition to LHC-B, a number of proponents from the former three collabo- 
rations decided not to continue. This unfortunately included P. Schlein and his UCLA 
colleagues who had pushed very strongly for the collider-mode idea with COBEX, and who 
had been a driving force up to that point. The LHCC decision sent a strong signal to the 
high-energy physics community: that CERN were prepared to give their strong support to 
one dedicated B experiment. This encouraged many physicists from institutions around 
the world to join the new collaboration over the subsequent years. T. Nakada, who was 


an instigator of the Gajet proposal, was elected LHC-B spokesperson. The 1995 Letter 
of Intent of LHC-B established the basis for the new detector design, which was refined 
in the following years, until its approval in 1998. By that time, the name had changed 
from LHC-B to simply LHCb. Fig. [1.3] shows the detector as it appeared in the Letter of 
Intent. 
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Figure 1.3: The LHC-B detector as it was proposed in 1995. All the basic components shown 
would be part of the final detector, albeit with many refinements and optimisations that will be 
described in the following Section. 


Soon afterwards, a competing B-physics experiment, BTeV [27], was proposed to run 
at the Fermilab Tevatron, incorporating a single magnet, a double arm spectrometer and 
a vertex trigger at the first level, in order to recover the reduced b production cross section. 
Following the decision to shut down the Tevatron, BTeV was not approved, however 
several of the experiment’s innovative ideas were carried through to the future LHCb 
experiment. 


2 The LHCb Detector 


The basic mechanism for heavy-quark production at the LHC is via gluon-gluon fusion. 
The angular distribution of cé or bb pairs is peaked at small angles with respect to the 
beam-line, with high correlation between the constituents of the pair. This allows the 
detection with good acceptance of the resulting hadrons in a rather limited solid angle. 
QCD calculations give cross-section values of cæ œ 1.5 mb and o œ 0.5 mb, respectively. 
The large event multiplicity requires a high granularity of the detector, together with 
minimal thickness in terms of radiation and interaction lengths to reduce secondary 
interactions. 

The LHC-B Lol presented a forward spectrometer with 400 mrad acceptance and 
a single large dipole magnet. The apparatus would be located inside the former DELPHI 
cavern at LEP with little modifications to the existing infrastructure but, for reasons of 
cost, the detector sacrificed half of the solid angle by only being a single arm spectrometer. 
The Lol design inherited important features from the three ancestor experiments and 
from the contemporary HERA-B, which had rate and radiation issues similar to those 
expected at the LHC. In contrast to the latter experiment, LHC-B had in addition a 
hadron calorimeter and a second (upstream) RICH with two radiators. Initially it was 
thought that an efficient tracking system in a harsh environment would require a large 
number of tracking stations, so, paralleling HERA-B, LHC-B had twelve tracking stations 
in the large-angle region. This number was reduced to ten in the Technical Proposal 
presented in 1998 [28], which by then had changed its name to LHCb. 

The disappointing performance of HERA-B was largely ascribed to the large amount of 
material in the detector, which prompted the LHCb collaboration to perform a thorough 
review of the apparatus, with the aim to reduce material without sacrificing performance. 
A Technical Design Report submitted in 2003 presented the LHCb “Reoptimized” Detec- 
tor (29). This is the basis on which the experiment was eventually built and is described 
in the following subsections. 


2.1 Overview 


The LHCb detector is a forward spectrometer, shown in Fig. and is installed 
at Intersection Point 8 of the LHC. A modification to the LHC optics, shifting the 
interaction point by about 11 m from the centre, allowed maximum use of the cavern 
space. This results in a detector length of approximately 20 m, and with maximum 
transverse dimensions about 6 x 5 m3. The angular acceptance ranges from approximately 
10 mrad to 300 mrad (Fig. in the horizontal magnetic-bending plane, and from 
10 mrad to 250 mrad in the vertical plane (Fig. [2.1p. With this geometry the detector is 
able to reconstruct approximately 20% of all bb pairs produced. 

To measure the momenta of charged particles, a dipole magnet producing a vertical 
magnetic field is used. It is a warm magnet providing an integrated field of 4 Tm, with 
saddle-shaped coils in a window-frame yoke, and with sloping poles in order to match 
the required detector acceptance. The design of the magnet allows for a level of fringe 
field inside the upstream Ring Imaging Cherenkov detector (RICH 1, see Sec. of less 
than 2 mT whilst providing a residual field in the regions between the upstream tracking 
stations. 

The tracking system consists of a silicon VErtex LOcator detector (VELO) (31), 


surrounding the interaction region, and four planar tracking stations, the TT tracker 
upstream of the dipole magnet and three tracking stations T1-T3 downstream of the 
magnet (321 |33). The T1-T3 stations consist of an Inner Tracker (IT), located at the 
centre of the stations and surrounding the beam-pipe, and an Outer Tracker (OT) for 
the outer regions. A minimum momentum of around 1.5 GeV/c is required for a track to 
reach the downstream stations [34]. 
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Figure 2.1: The LHCb detector: side view 
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Figure 2.2: The LHCb detector: top view 


Particle identification (PID) is a fundamental to the goals of the LHCb experiment 
by separating pions, kaons and protons produced in heavy-flavour decays. Accurate 
reconstruction of electrons and muons is crucial for flavour tagging. All aspects of PID 
are accomplished by a set of specialized detectors. 

The PID system is based around two RICH detectors, designed to cover almost the 
full momentum range of tracks in LHCb. The upstream detector, RICH 1, covers the 
low momentum region from œ 2 to 60 GeV/c using aerogel (in Run 1 only) and Cy4F jo 
radiators, whilst the downstream detector, RICH 2, covers the high momentum range 
from 15 GeV/c up to and beyond 100 GeV/c using a CF, radiator. 

Two calorimeters, one electromagnetic (ECAL) and the other hadronic (HCAL), 
supplemented by a Preshower Detector (SPD/PS) [35], provide identification of electrons, 
photons and hadrons and a measurement of their energy. This measurement is used at 
the trigger level to select candidates on the basis of their transverse energy. Muons play a 
crucial role in many of LHCb’s measurements because of the cleanliness of the signature. 
Their identification is achieved by five muon stations (M1 — M5), interspersed with iron 
filters. The muon system also supplies measurements of muon transverse momenta for the 
trigger. 


2.2 The Vertex Locator 


The role of the VELO is to measure the impact parameters of all tracks relative to the 
primary vertex (PV), to reconstruct the production points and decay vertices of hadrons 
containing b- and c-quarks and to allow precision measurements of their mean lifetimes. 
The subdetector accepts particles with pseudorapidities in the range 1.6 < 7 < 4.9 and 
which have PVs within |z| < 10.6 cm from the nominal collision point along the beam 
direction. The VELO is split into two halves surrounding the beam-pipe, each containing 
21 modules. Each module is then made up of two silicon half discs of 300 um thickness, 
one with strips in the radial, r coordinate, the other in the polar, @ coordinate. This 
cylindrical geometry allows a fast track- and vertex-reconstruction to be made at the 
second stage of trigger. The strip segmentation is such to limit the highest occupancy of 
the strips to less than 1.1 %. 

The VELO is positioned, with an accuracy better than 4 um, at the closest distance 
possible from the beam, about 7 mm during data taking, The sensors operate within 
a so-called Roman pot configuration, located inside a secondary vacuum of less than 
2.1077 mbar pressure, separated from the primary LHC vacuum. The sensors are 
retracted during beam injection and are then quickly moved in for physics operation when 
the LHC beams are stable. 

The vessel containing the silicon discs and the front-end electronics (RF-box) has 
aluminium walls of 300 wm thickness to minimize multiple scattering. The average 
material budget of the detector for tracks in the LHCb acceptance is 0.22 Xo. In order to 
minimize radiation damage and to dissipate the produced heat, a cooling system keeps 
the temperature range between -10 to 0°C. 

Fig. [2.3] summarizes the VELO performance in terms of impact parameter and decay 
time resolution (31). 
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Figure 2.3: Left: Projected impact parameter resolution as a function of 1/pr. Right: Decay 
time resolution (points) as a function of momentum for B}? —> J/@ decays. The superimposed 
histogram shows the distribution of momentum for the decay. 


2.3 The TT and Downstream Tracking System 


Following the VELO, the tracking system is composed of the TT station, located between 
RICH 1 and the magnet, and three stations (T1,T2,T3) downstream of the magnet. The 
TT is composed of four stations grouped in pairs, called TTa and TTb, spaced by 30 cm. 
Each station consists of silicon microstrip planar modules covering a rectangular area of 
150 cm x 130 cm (width times height), covering the LHCb acceptance of 300 mrad in 
the horizontal plane and 250 mrad in the vertical. The strips of the first and the fourth 
stations are vertical and measure the bending x coordinate, whilst the second and third 
planes have stereo angles of +5°, respectively. 

Tracking stations T1 - T3 each consist of an inner part (IT) surrounding the beam 
pipe, and an outer part (OT) beyond. Each IT station consists of four overlapping silicon 
layers, two rotated by a stereo angle of +5 deg and two aligned to the y (vertical) axis. 
Each layer is made up of four independent modules placed around the beam pipe, covering 
about a 120 x 40 cm? area, as shown in Fig. 


Figure 2.4: The LHCb inner detector, (left) a vertically-aligned layer, (right) a stereo layer. 
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The spatial resolutions of both the TT and IT are approximately 50 um per hit, with 
strip pitches of about 200 um. The hit occupancies vary between 1.9% for the inner 
sectors to 0.2% for the outermost modules. To minimize radiation damage, the sensors 
operate at 5°C temperature. 

The OT is a drift detector consisting of straw tubes with internal diameters of 4.9 
mm, each filled with an Ar/CO, gas mixture in a ratio 70 — 30%. The straws provide 
a 35 ns maximum drift time and 205 um spatial resolution with 17% maximum straw 
occupancy. Each of the three stations is made of four modules, shown schematically in 
Fig. [2.5] A picture of the assembled OD is shown in Fig. [2.6] In the first and third modules 
the ee tubes are aligned to the vertical axis while the third and fourth modules have 
stereo angles of +5 deg. The total active area is about 5.97 x 4.85 m?, covering the full 
LHCb acceptance. 

The overall tracking efficiency for “long” tracks (i.e. those tracks measured in all the 
tracking detectors including the VELO) is greater than 96 % for 5 < p < 200 GeV/c. The 
momentum resolution dp/p is 0.5% at low momentum, increasing to 1.1% at 240 GeV/c. 
The mass resolution is 14.3 MeV/c? for the J/) resonance. 


eh flea MTT yi iii 1) 


e 


i 


Figure 2.5: Schematic views of the LHCb outer detector. 


2.4 The RICH system 


The role of the RICH system is to provide m/K/p discrimination for LHCb, which is 
essential for most CP-violation studies, background rejection and flavour tagging. The 
momentum range which contains 90% of kaons, pions and protons from B meson decay 
is between 2 and 150 GeV/c and, to achieve this separation, two Cherenkov detectors, 
RICH 1 and RICH 2, are employed. 

The RICH 1 detector differentiates particles with low and intermediate momenta, from 
1 to ~60 GeV/c. It is located close to the interaction region, upstream of the magnet, and 
covers the acceptance from +25 mrad to +300 mrad (horizontal plane) and to +250 mrad 
(vertical plane). RICH 1 initially contained two different radiator materials: an aerogel 
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Figure 2.6: The LHCb Outer Detector in place on the beam line. Well visible is the beam pipe. 


layer 5 cm thick with refractive index n = 1.03 and a C4F 9 gas layer of length 85 cm with 
refractive index n = 1.0014. Aerogel has the power to provide r/K discrimination from 
about 1 up to 10 GeV/c, however it was removed for Run 2 due to occupancy problems. 
The C,Fj9 radiator extends the positive 7/K identification from about 10 GeV/c to 
60 GeV/c, however r/K discrimination below 10 GeV/c is still possible by operating the 
RICH in kaon veto mode. 

RICH 2 has a smaller angular acceptance of +15 mrad to £120 mrad (horizontal plane) 
and to +100 mrad (vertical plane) and covers the region where high momentum particles 
are most abundant. It is located downstream of the magnet, between T3 and the first 
muon station M1. RICH 2 uses a CF; gas radiator with refractive index n = 1.00046. 

In both detectors, Cherenkov photons are detected by a combined system of plane 
and spherical mirrors to focus photons onto a pair of photo-detector planes, where Hybrid 
Photon Detectors are employed to detect the Cherenkov rings. The photo-detectors are 
located outside the detector acceptance, in regions of low magnetic field and relatively 
low radiation. 

The Cherenkov angles for the three different RICH radiators and for different particles 
as a function of momentum are shown in Fig. |2.7| (left). A measurement of RICH 
performance in LHCb data is shown for the two gaseous radiators In Fig. (right) [36]. 


2.5 Calorimeters 


The calorimeter system identifies hadrons, electrons and photons, and also measures 
their energies and positions for the Level-0 trigger. The system is composed of four 
sub-detectors: the scintillator Pad Detector (SPD), the PreShower detector (PS), the 
Electromagnetic Calorimeter (ECAL), and the Hadron Calorimeter (HCAL). The SPD 
and PS are located just upstream of the ECAL. The ECAL, PS and SPD are segmented 
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Figure 2.7: (Left) Values of Cherenkov angle as a function of momentum for different particles 
for the three RICH radiators of refractive index n. (Right) Measured Cherenkov angles in LHCb 


data [36]. 


into three sections in the xy plane, with active pads growing from the inner to the outer 
regions. The HCAL is similarly divided in two sections. The corresponding granularities 
are outlined in Fig. 


Figure 2.8: The (left) LHCb electromagnetic calorimeter and (right) hadronic calorimeter. The 
bottom left regions indicate the areas occupied by the beam-pipe. 


The SPD and PS are used at the trigger level and offline, in association with the 
ECAL, to indicate the presence of electrons, photons and neutral pions. The detectors 
have two plastic scintillator layers separated by a 15 mm thick lead plate where electrons 
and photons can radiate; the downstream scintillator then samples the radiated energy. 
The light from the scintillators is sent to photomultipliers by wavelength-shifter (WLS) 
optical fibers. 

The ECAL employs the Shashlik technology, where independent modules, constructed 
from scintillating tiles and lead plates, are alternated (see Fig. [2.9). The ECAL has 
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66 layers of such modules consisting of 2 mm of lead followed by 4 mm of scintillator 

material. The ECAL also uses WLS optical fibers to guide the light from the detector 

to photomultipliers, placed on the back face of each module. The energy resolution 
achieved is 

o(E) 10% 

E VE 


where F is the electron energy expressed in GeV. 

The HCAL is also a sampling calorimeter. It is constituted of iron absorber with 
scintillating tiles as the active material. The innovative feature of this sampling structure 
is the orientation of the scintillating material: the tiles run parallel to the beam axis. In 
the lateral direction, tiles are spaced with 1 cm iron, while longitudinally the length of 
the tiles and iron spacers correspond to the hadron interaction length A; ~ 20 cm in steel. 
Light is collected by WLS optical fibres running along the detector towards the back side 
where the photomultiplier tubes are located (see Fig. [2.9). 


91% (2.1) 


Z4 Wavelength-Shifting Fibres 


Scintilator Tile —— 


Figure 2.9: A schematic showing the layout and segmentation of the LHCb ECAL and HCAL 
calorimeters. 


The HCAL is used to measure the hadronic shower transverse energy for the Level-0 
trigger and to improve the high momentum electron/hadron separation. The energy 
resolution achieved is . 69 5)% 

CE _ OG 15 0) (2.2) 
E VE 


where F is the hadron energy in GeV. 


2.6 The Muon System 


The Muon System consists of five stations, M1-M5, of rectangular shape. The complete 
system is made up by 1368 Multi Wire Proportional Chambers supplemented by 12 Triple 
GEM Chambers in the inner region of the first station, to cope with the very high particle 
rate. The chambers employ a variety of readouts, optimized for a precise pr measurement 
for the trigger. The complete system has an acceptance in the bending plane from 20 mrad 
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to 306 mrad, and in the non-bending plane from 16 mrad to 258 mrad. This results in a 
total acceptance of about 20% for muons from semileptonic inclusive b decays. 

The M1 station is located in front of the calorimeters and is used in order to improve 
the p; measurement for the trigger. The geometry of the five stations is projective; all the 
transverse dimensions scale as the distance from the interaction point. Stations M2-M5 
are placed downstream of the calorimeters and are interleaved with 80 cm thick iron 
absorbers. The total absorption thickness, calorimeters included, is about 20 interaction 
lengths. In this way the minimum momentum for muons crossing the five stations is about 
6 GeV/c. 

Each muon station is designed to achieve an efficiency above 99% in a 20 ns time 
window with a noise rate below 1 kHz per physical channel, as described in [37]. To reach 
such an efficiency, four chamber layers per station are used in M2-M5 (two layers in M1). 
The time resolution is achieved by a fast gas mixture Ar/CO2/CF, in the ratio 40:55:5. A 
ratio 45:15:40 is employed in the Triple GEM chambers. 


2.7 The trigger 


Even with the relatively large bb cross-section at LHC energies, only approximately 1% 
of visible pp interactions result in a bb event. Moreover, only about 15% of those events 
will produce at least one b-hadron with all decay products passing within the acceptance 
of the spectrometer. The branching fractions of decays used to study CP violation are 
typically less than 1073. Further reductions are unavoidable in the offline selection, where 
stringent cuts must be applied to enhance signal over background. Therefore the purpose 
of the LHCb trigger is to achieve the highest efficiency for the events later selected in the 
offline analysis while rejecting drastically most of the uninteresting background events. 
To achieve this goal, the trigger uses information from all LHCb sub-detectors. 

The trigger is organised in two different levels: the Level-0 (LO) trigger based on 
custom electronic boards, and the High-Level Trigger (HLT), implemented in a computer 
farm. Level-O uses the information from the calorimeter and muon systems, performing 
a selection in order to reduce the event rate from 40 MHz to below 1 MHz, which is 
the maximum frequency allowed to read out the entire detector. The HLT is a software 
application running on a processor farm that further reduces the rate of events in the kHz 
range for storage (see Fig. |2.10). 

The HLT has significantly evolved over time from the original design in the LHCb 
Technical Proposal (TP) in 1998, to the trigger design in the Technical Design Report 
(TDR) in 2003, to the Run 1 (2010-2012) actual implementation and finally to the 
additional features introduced during Run 2 (2015-2018) (41). In the TP it was assumed 
that a first HLT trigger level (L1) would reduce the 1 MHz input rate to a 40 kHz output 
rate with a variable latency of less than 256 us, using coarse information from the vertex 
detector to reconstruct vertices and tracks with no momentum information (the VELO r-¢ 
geometry was designed for this purpose). A second HLT trigger level (L2) was fashioned 
to extrapolate VELO tracks into the magnetic field to the tracking stations downstream of 
the magnet and reduce the output rate to 5 kHz with an average latency of 10 ms. Finally 
a third level (L3) would implement the full event reconstruction and a set of exclusive 
selections to bring down the rate to 200 Hz. 

By the time of the trigger TDR in 2003, it became clear that the LHC was not going to 
start before the end of the decade when much more powerful processing units would become 
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available. In addition a series of test-beam and detailed simulation studies convinced the 
collaboration of the need to have momentum information at the first stage of the HLT. 
Therefore a new tracking station just upstream of the magnet was introduced (the TT 
station). In addition, a shield which had been protecting RICH1 from stray magnetic fields 
was removed to allow for a rough estimation for the momentum of tracks reconstructed 
between the VELO and TT stations. The software trigger then had two levels: Level-1 
able to reduce the output rate to 40 kHz using LO, VELO and TT information with an 
average latency of 1 ms, and HLT to reduce the output rate to 200 Hz with a combination 
of inclusive and exclusive selections. Between the time of the trigger TDR and the first 
physics run (Run 1), the interest in having a more performing HLT for charm physics 
(c@ with a factor 20 larger production cross section than bb) and a much more robust 
system, convinced the collaboration to push for much more inclusive selections in the final 
trigger stage and a much larger trigger output rate (3-5 kHz). This implied a complete 
redefinition of the offline data processing model. 

Furthermore, it had been assumed that the LHC would operate with a 25 ns bunch 
separation, limiting the number of overlapping events to a mean number of u ~ 0.4 per 
bunch crossing at a luminosity of 2 x 10°? cm~?s~!. When, from 2011, a separation of 
50 ns was adopted for early LHC operation, the experiment decided to run at u ~ 1.4 
to compensate for the lower number of bunches. Therefore the HLT had to adapt to 
running conditions rather different than first assumed. This was made possible by the 
highly flexible design of the HLT. 

After the success of the LHCb trigger performance in Run 1, the good understanding 
of the trigger reconstruction allowed the introduction of the ” real-time analysis” concept 
during Run 2. After the first HLT trigger level (HLT 1), events are buffered to disk storage 
in the online system. This is done for two purposes, firstly events can be processed 
further during inter-fill periods, and secondly the detector can be calibrated and aligned 
run-by-run before the HLT2 stage. Once the detector is calibrated and aligned, events are 
passed to HLT2, where a full event reconstruction of “offline quality” is performed. This 
allows for a wide range of inclusive and exclusive final states to trigger and obviates the 
need for further offline processing. In addition, new techniques to reduce the amount of 
information saved per event allowed to increase significantly the output rate to 10-15 
kHz, as in Fig. while the output of HLT1 could be increased to O(110 kHz). The 
decrease in requests for offline reconstruction also helped to mitigate the pressure on the 
offline computing model. 


2.7.1 Level-0 hardware trigger 


The LO trigger is divided into three independent components: the LO-Calorimeter trigger, 
the LO-Muon trigger and the LO-PileUp trigger. The latter is used to reject multiple 
visible interactions in a bunch crossing by means of the “ad hoc” Pile-Up System detector 
housed in the VELO. The first two components are briefly described below. 


The LO-Calorimeter part of the trigger obtains informations from the SPD, PS, ECAL 
and HCAL subdetectors and computes the transverse energy deposited by incident particles: 
Er = Eocos@, where Ep is the energy of the particle and @ is the polar angle given by 
the cell hit in the detector. Together with energy information, the total number of hits in 
the SPD (SPD multiplicity) is also determined in order to veto large multiplicity events 
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Figure 2.10: Trigger overview in Run 1 (left) and Run 2 (right) 


that would take too large a fraction of the available processing time in the HLT. From 
the calorimeter information, three types of candidates are built and selected according to 
specific Er criteria: i) Hadron candidate (LOHadron); ii) Photon candidate (LOPhoton); 
and iii) Electron candidate (L0Electron). 


The LO-Muon part of the trigger requires a muon candidate to have a hit in all five 
muon stations. The LO muon processor boards select the two highest pr muon tracks in 
each quadrant of the muon system with a maximum of eight candidates. The trigger sets 
a single threshold either on the largest muon pr (LO muon trigger) or on the product of 
the largest and the 2nd largest (LO dimuon trigger). Events with SPD multiplicity > 600 
are excluded in the LO muon trigger in order to minimize the track multiplicity. This 
limit is raised to 900 in the LO dimuon trigger at the expense of a small increase in rate. 

The total output rate of the LO trigger is limited to 1 MHz, which is the maximum 
rate accepted by the HLT1. Such an output rate consists of about 400 kHz of muon 
triggers, about 450 kHz of hadron triggers and about 150 kHz of electron and photon 
triggers (the individual triggers have an overlap of about 10%). 


2.7.2 High Level Trigger 


Data from LO are sent to the Event Filter computer Farm (EFF) which runs the HLT 
algorithms. The HLT is a software application whose 29500 instances run on the EFF. 
Each instance is made up of independently operating trigger lines; each line consists of 
selection parameters for a specific class of events. 

The HLT is divided into two stages. The first stage (HLT1) processes the full LO rate 
and uses partial event reconstruction to reduce the rate to about 110 kHz. The second 
stage (HLT2) reduces the rate to about 12.5 kHz, performimg a more complete event 
reconstruction [41]. 
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HLT1 reconstructs the trajectories of charged particles traversing the full LHCb 
tracking system which have a pr larger than 500 MeV. The hits in the VELO are 
combined to form straight-line tracks loosely pointing towards the beam line. Next, at 
least three hits in the TT are required in a small region around a straight-line extrapolation 
from the VELO. The TT is located in the fringe field of the LHCb dipole magnet, which 
allows the momentum to be determined with a relative resolution of about 20%, and this 
estimate is used to reject low pp tracks. Tracks are then extrapolated to the T-stations 
downstream of the magnet. The search window in the IT and OT is defined by the 
maximum possible deflection of charged particles with pr larger than 500 MeV. The search 
is also restricted to one side of the straight-line extrapolation by the charge estimate of 
the track. Subsequently, all tracks are fitted with a Kalman filter to obtain the optimal 
parameter estimate using a simplified geometry description of the LHCb detector. The 
set of fitted VELO tracks is re-used to determine the positions of the PVs. 

Tight timing constraints in HLT1 mean that most particle-identification algorithms 
cannot be executed. The exception is muon identification due to its clean signature. Hits 
in the muon stations are searched for in momentum-dependent regions of interest around 
the track extrapolations. Tracks with p < 3 GeV cannot be identified as muons, as they 
would not be able to reach the muon detectors. 

HLT1 has two inclusive trigger lines which select events containing a particle whose 
decay vertex is displaced from the PV: a line which selects a single displaced track with 
high pr, and a line which selects a displaced two-track vertex with high pr. Both lines 
start by selecting good quality tracks that are inconsistent with originating from the PV. 
The single-track trigger then selects events based on a hyperbolic requirement in the 2D 
plane of the track displacement and pr. The two-track displaced vertex trigger selects 
events based on a multivariate discriminant whose input variables are the vertex-fit quality, 
the vertex displacement, the scalar sum of the pr of the two tracks and the displacement 
of the tracks making up the vertex. The two-track line is more efficient at low pr, whereas 
the single track line performs better at high pr, such that in combination they provide 
high efficiency over the full pr range. 

The HLT1 muon lines select muonic decays of b and c hadrons, as well muons originating 
from decays of W and Z bosons. There are four main lines: one line that selects a single 
displaced muon with high pr, a second single muon line that selects very high pr muons 
without displacement for electroweak physics, a third line that selects a dimuon pair 
compatible with originating from a decay of a charmonium or bottonium resonance or 
from Drell-Yan production, and a fourth line that selects displaced dimuons with no 
requirement on the dimuon mass. During Run 2, typically about 80 kHz were allocated 
to the inclusive HLT1 lines, while about 20 kHz to the muon lines. The rest of the HLT1 
output is dedicated to special low multiplicity triggers and calibration trigger lines. 


HLT2 can perform the full event reconstruction since the output of HLT1 is buffered. 
The full event reconstruction consists of three major steps: the track reconstruction of 
charged particles, the reconstruction of neutral particles and particle identification. The 
HLT2 track reconstruction exploits the full information from the tracking sub-detectors, 
performing additional steps of the pattern recognition which are not possible in HLT1. 
Tracks with a pr larger than 80 MeV are reconstructed in HLT2, without the requirement 
to have hits in the TT station. This is to avoid inefficiencies due to the TT acceptance, 
which is crucial for part of the charm and kaon physics programme. In addition, tracks 
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produced by long-lived resonances that decay outside the VELO are reconstructed using 
T-station segments that are extrapolated backwards through the magnetic field and 
combined with hits in the TT. Similarly, the most precise neutral cluster reconstruction 
algorithms are executed. Finally, in addition to the muon identification available in HLT1, 
HLT2 exploits the full particle identification from the RICH detectors and calorimeter 
system. 

The HLT2 inclusive b-hadron trigger lines look for a two-, three-, or four-track vertex 
with sizeable pr, significant displacement from the PV, and a topology compatible with 
the decay of a b-hadron, using a multivariate discriminant. Whenever one or more tracks 
are identified as muons, the requirements on the discriminant are relaxed to increase the 
efficiency. As in the case of HLT1, several muon lines are used to select muonic decays of 
b and c hadrons and of W and Z bosons. However in HLT2, the muon reconstruction is 
identical to the offline procedure, having access to exactly the same information. During 
Run 2, typically about 3 kHz of the trigger rate is from the inclusive b-hadron trigger 
while the muon lines take about 1 kHz. A large fraction of the trigger bandwidth (2-4 kHz) 
is allocated to exclusive selection of charm decays, where a reduced amount of information 
is saved per event. The rest of the trigger bandwidth is due to other special triggers and 
calibration trigger lines. 
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3 LHCb contributions to CKM measurements and 
CP violation 


The violation of the combined operation of charge conjugation and parity, CP, was first 
observed in 1964 in decays of neutral kaons [43]. The BaBar and Belle B Factory 
experiments and the CDF experiment established CP violation in the decays of neutral 
B? mesons. LHCb now extends measurements to much greater precision, and also probes 
the B, system, which is vital to explore the full range of CP violation measurements. 

In the Standard Model, the Cabibbo-Kobayashi-Maskawa (CKM) unitary matrix ajen], 
Vcxm, describes the electroweak coupling strength V;; of the W boson to quarks 7 and j: 


Vaa Vus Vab 
Vekm = | Vea Ves Væ |- (3.1) 
Via Vis Vio 


CP is violated in the Standard Model if any element of the CKM matrix is complex. 
The parametrisation of the CKM matrix due to Wolfenstein is given by 


a $ AX? (p — in) 
Vorm = | -A+ 34A- 2o + in] 1 — 3? - 341 +44?) AN 
Ad[L — (1— 42)(p + in)] —AX? + 4AM [1 — 2(p + in)] 1 — 442M 2 
3.2 


for the four Standard Model parameters (A, A, p,7)). The expansion parameter, À, equal to 
the sine of the Cabibbo angle, has a value |V,,,| = 0.22 [49], and in Equ.[3.2] the expansion 
is given for terms up to order \°. 

The unitarity of the CKM matrix leads to six orthogonality conditions between any 
pair of columns or any pairs of rows of the matrix. The orthogonality means the six 
conditions can be represented as six triangles in the complex plane. The interesting 
relations for CP violation are those given by: 


Vaa Vi + VeaVn + VaV =0 : the unitarity triangle, (3.3) 
us Vay + Ves Va + VisVa = 0 : the B, triangle, and (3.4) 
Vea Vga + VerViin + Ves Vis =O : the charm triangle. (3.5) 


The unitarity triangle has sides with lengths that are the same order in A, namely O(A°), 
which implies large CP asymmetries in B° and B* decays. The B, triangle has two sides 
of O(A?) and the third of O(\*). Hence CP violation in B, mixing is significantly smaller 
than in the B° system. Moreover, the charm triangle has two sides of O(A) and the third 
of O(A°), hence CP violation in the charm system is expected to be extremely small. Note 
that all three triangles have equal area [50]. 

To study CP violation, the B—physics experiments measure the complex phases of 
the CKM elements and measure the lengths of the sides of the triangles to check for a 
self-consistent picture. CP violation is predicted in many (often very rare) B hadron 
decays, hence LHCb utilises large samples of B, B,, Be mesons and B—baryons. New 
physics can be discovered and studied when new particles appear in, for example, virtual 
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Figure 3.1: The knowledge of the unitarity triangle as of autumn 2009 [52]. 


loop processes of rare B decays, leading to observable deviations from Standard Model 
expectations, both in branching ratios and CP observables. Hence the LHCb strategy is 
to determine with high precision the CKM elements and to compare measurements of the 
same parameters, especially those where one is sensitive to new physics and the other to 
Standard Model processes. 


3.1 The status of the unitarity triangle before LHCb 


The first generation B Factory experiments to study CP violation in the B—system, 
BaBar and Belle, made huge in-roads into testing the Standard Model description of CP 
violation; the status was summarised extensively at the Beauty 2009 Conference (51). 
CDF and DO extended these studies at the Tevatron, and make first explorations in the 
B, sector. Fig. [3.1]shows the status of the unitarity triangle measurements compiled by 
the CKM-Fitter Group in 2009, when the B Factories had been running for around 
ten years. Here graphical results are displayed in the p — 7 plane and the best fit to the 
apex of the triangle (Equ. to the 95% confidence level is shown. The fit to the CKM 
parametrisations include measurements of the sides of the triangle through measurements 
of the CKM elements and the angles, information from rare K and B meson decays, and 
B? — Bo mixing. 

Before 2009, when the LHC turned on, the B Factory experiments, CDF and, to a 
lesser degree DO, measured the parameters of the unitarity triangle with varying degrees 
of precision: 


e The quantity sin 23 was measured in all channels, including the “gold plated” channel 
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B® + J/y K$, to a precision of around ~0.03; 


e The sides |Via/Vis| and |Vi/Vs| were known from B? — B? mixing and from b > u 
decays, respectively each to ~10%, but limited by theory. The B, mixing phase (@,) 
was unmeasured; 


e The angle a was measured in the channels B > a7, pr and pp with a statistical 
precision of ~ 5°; 


e There was a statistics-limited measurement in B + DK modes of the angle y to 
around 20 — 25°. A measurement of y from B, modes such as B? + Df K~ had 
been completely unexplored. 


e The parameter €g, measured in kaon decays, provided a very loose constraint on 
the triangle vertex; 


e The B, sector had been largely unexplored by the first generation experiments, as 
had b—baryons. Running at the Y(45), the B Factories produced predominantly 
Bu a meson pairs, however Belle record a significant sample of T(55') data, allowing 
some interesting measurements of B, pairs. At the LHC, B°, B*, B., Be and 
b—baryons are produced in approximately in the ratios ~ 40 : 40 : 10: 0.1 : 10%; 


e The B Factories were statistics limited for very rare processes with branching ratios 
< 1x 10~°, such as b > s flavour-changing neutral current (FCNC) transitions, 
e.g. b > sy and b > sll. Super-rare transitions such as B(s, > ut were also 
unobserved. 


In contrast, Fig. shows the status of the unitarity triangle measurements today 


Figure 3.2: The current knowledge of the unitarity triangle (as of Summer 2019) [52]. 
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Figure 3.3: The box diagrams representing mixing between B? and By mesons. 


3.2 Heavy quark mixing measurements 


Since flavour is not conserved in the weak interaction, mixing between B° and B° mesons 
(where q = d or s) is possible via the box diagrams shown in Fig. |3.3} The probability 
for finding a B, (or a B,), given the initial state was a B, (or a B,) at time At after 
production, is given by: 


Pmiz( At) = exp Žž (1 — cos(Am,At)) . (3.6) 

Here Am, is the mass difference my — mz, where myz are the masses of the heavy and 

light mass eigenstates and Tpo is the B lifetime. At the LHC, the two neutral B mesons 
produced can oscillate independently at any time after production. 

Any CP measurement from a time-dependent analysis of neutral B decays needs the 

determination of the B flavour (b or b) at production. This requires b—quark “tagging”, 


and several algorithms have been developed by LHCb involving the combination of 
so-called opposite side and same side taggers (54). 


e B® mixing: B° oscillations are measured in the channel B° > D~y*v,,X and the 
charge conjugate modes. The oscillation measurements rely on tagging, and are 
shown for 3 fb~! of data in Fig. (left). The LHCb Am, measurement, which is 
the world’s best, is Am4 = (505.0 + 2.1 + 1.0) ns“ [55]. 


e B, mixing: B? — B? oscillations were first observed by CDF [e]. B, oscillations 
are measured at LHCb in the mode B}? + D7r* and its charge conjugate state. The 
oscillations are shown for 1 fb~! of data in Fig.|3.4] (right). The plot shows the proper- 
time distribution of B? + D;7* candidates, in five different B? decay channels, that 
have been flavour-tagged as not having oscillated. The LHCb Am, measurement, 
which is again the world’s best, is Am, = (17.768 + 0.023 + 0.006) ps7! [57]. 


e Charm meson mixing: 


D? -D° mixing was established in a single experiment in 2013 by LHCb in D® > 
K*n~ decays, although it was confirmed earlier by combining several B Factory 
results. The D? mixing parameters were measured in LHCb in a decay-time- 
N = : . í : 
dependent fit to the ratio oe The time-dependent fit is shown in Fig!3.4 
DO4K~at+ 


(lower) for 1 fb~! of data. The no-mixing scenario is excluded at 9.lo |58]. 
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e Tagged mixed 
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Figure 3.4: Proper-time distribution of (left) B° — B’ oscillations in B? > D- utv, X decays, 
(right) B? — B oscillations in B? + D7 r* decays, and (lower) D? — D’ oscillations in D° > 
K*n~ decays. All distributions rely on flavour tagging, and the curves correspond to the fitted 
oscillations. 


3.3 LHCb measurements of the unitarity triangle parameters 


The LHCb experiment performs a high-statistics study of CP violation with unprecedented 
precision in many different and complimentary channels, providing a sensitive test of the 
Standard Model and physics beyond it. 


3.3.1 Measurements of the CKM angle 6 


The time-dependent decay asymmetry of the channel B? > J/y K} allows a measurement 
of the angle 6. This is known as the “golden” decay mode because the channel is virtually 
free of penguin pollution (which enters with the same overall phase), resulting in very 
small theoretical uncertainty, of order 1% (59). CP violation in this channel occurs in the 
interference between mixing and decay, where the mixing process introduces a relative 
CP-violating weak phase of 26. Experimentally the CP asymmetry is measured from the 
ratio of the numbers of B and B mesons, Ng, y and Np-,r, decaying into final state f: 


Acp = Seat Neos = sin (Amat) sin (28) — cos (Amat) cos (28). (3.7) 
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The LHCb measured and fitted asymmetries for the J/w (1S) and (2S) states are 
shown in Fig. for 3 fb! of data at 7 and 8 TeV (60). These measurements are 
cos (Amat) = —0.017 + 0.029 and sin (Amat) = 0.760 + 0.034, where an observation of 
a direct CP-violation contribution proportional to cos (Amat) would be an indication 
of new physics [59]. The LHCb measurement is now competitive with BaBar and Belle 
measurements; the current world average of sin28 = 0.695 + 0.019 is dominated 
by LHCb together with the B Factory measurements in the complementary channels 
B® + J/WK3 and B? + J/wWK{. The measurement by LHCb of sin2@ in gluonic 
penguins will further contribute to this study. 
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Figure 3.5: The LHCb measured and fitted asymmetry in B? + J/WK decays for (left) the 
J/wW(1S) and (right) the 7(2S). 


3.3.2 Measurements of the CKM angle a 


The primary method at LHCb for the measurement of œ is through an amplitude analysis 
via the B — pr decay modes [62], however these channels are difficult at LHCb due 
to the need to efficiently reconstruct 7°s. Penguin pollution is present and must be 
constrained, with the additional application of isospin symmetry. The precision on Q at 
LHCb is expected to be dominated by systematic uncertainties, and any measurement 
is not expected to improve on a combination of the B Factory measurements, a = 


(86.4435) [52]. 


3.3.3 Measurements of the CKM angle y 


A precise measurement of the angle y is key to understanding the closure (or otherwise) 
of the unitarity triangle. Constraints on the unitarity-triangle apex largely come from 
loop decay measurements which are very sensitive to the presence of new physics. ~y is the 
only angle accessible at tree level and hence forms a SM benchmark to which the loop 
measurements can be compared (assuming no significant new physics in tree decays). The 
y measurement also relies on theoretical input which is very well understood (63}/64). 

Determination of y from a combined fit to all measured parameters of the unitarity 
triangle currently gives a value y = (65.8447) [52]. Conversely the measurement of y 
alone from a combination of all direct measurements from tree decays gives y = (72.1727) $ 
Hence reaching degree-level precision from direct y measurements is crucial. 


25 


LHCb makes measurements of y by a variety of methods, where complementary is 
vital. Examples of the most sensitive LHCb measurements are outlined below. 


e y in the “time integrated” B+ — D°K+ modes 


The measurement of y is made in direct CP-violation via B* + D° K+ by three 
different methods: the GLW method (decay into a CP eigenstate) (65}(66}, the ADS 
method (decay into a flavour-specific mode) [67], and the GGSZ method (Dalitz 
analysis) fes]. These all access y through interference between the B+ — D° K+ and 
B+ > D K+ decay paths, where the D° and D decay to the same final state. When 
using these methods, the decay modes are self-tagging. In addition time-dependent 
analyses are not necessary. For the ALD and GLW modes, the charge-conjugate 
event yields are simply counted to determine the CP asymmetries (i.e. effectively a 
“counting experiment” ). 


Figure|3.6|shows an example of a CP asymmetry in the B+ + D°K* ADS mode [69]. 
Here the D°? is produced in a Cabibbo favoured mode (Væ) but decays via a 
suppressed mode (V.q) into Ktz~. This interferes with the D charge-conjugate 
state which is produced in a suppressed mode (V,,,) but decays to the same final 
state Kt via a favoured mode (V.,). The branching fraction for the favoured B 
decay is only ~ 1074, so these measurements require high statistics. The asymmetry 
observed in Fig. [3.6] has a magnitude of around 40% and has a significance of 7o. 
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Figure 3.6: Example of a CP asymmetry in the ADS mode B+ — D°K*, where the D° and D? 
both decay into K+*r~. The event yield is shown as a function of DK mass. 


A specific example of an LHCb analysis using the GGSZ Dalitz method is in the 
decay B+ > D°K* where D? + K$n+m- or D° > K}K+K- [70]. There is a rich 
Dalitz plot structure with the presence of large interference effects. The Dalitz space 
is divided up into symmetric bins, chosen to optimise sensitivity. An amplitude 
analysis can then be used to extract y. 


In all B* > D°K* modes and decays listed above, y can also be extracted from the 
corresponding B? + D°K** modes, albeit with reduced y sensitivity. In addition 
B? + D“ K“ GGSZ modes are also included in the global fit to extract the 7 
average value. 


e y from the “time-dependent” B? — D7 K+ mode 


The channel B? + D7 Kt, and its charge conjugate states, provide a theoretically 
clean measurement of the angle (y + ¢,) where @, is the (small valued) B, mixing 
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phase, with no significant penguin contribution expected (71). Here both B®, and 
via the mixing diagram B , can decay to the same final state D7 K*, resulting 
in interference which is sensitive to y. The same is true for decay into the charge 
conjugate state DIK. Hence four time-dependent decay rates are measured: 
B? + D>K+, B? + D}K-, B? > D7K+* and B? > Dt K~. The method is 
then to fit two asymmetries of the form 


N31 = Neat 


Acp = : 
a Nos + Neos 


(3.8) 


These measurements yield values for the strong phase difference dgcp between the 
amplitudes B > f and B —> f, the amplitude ratio, and (y + ¢,). 

The current measurement by LHCb in 1 fb~! of data yields a value y = ( 1 15°35) i ' 
This complements the measurements in B+ + D° K+, although with less statistical 
precision. 


e The y combination 


The LHCb measurement of y averaged over all the above methods, which includes 
all B°, B+ and B, modes, is y = (74.0855) [73]. This measurement dominates the 
current world average. The confidence limits as a function of y for the combination 
is shown in Fig. [8.7]for the various measurement channels. The agreement between 
B, and B* initial states is currently at the 20 level. 


WN B? decays 
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HE Combination 
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Figure 3.7: Confidence limits as a function of y in the LHCb combination for the various 
measurement channels. 


3.3.4 The sides of the triangle 


e The side opposite to 3 


Currently the closure test of the unitarity triangle is limited mainly by the side 
opposite to 8 which has a length proportional to |V,.|/|V.»| in the Standard Model. 
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This limitation is a consequence of tension between B Factory inclusive and exclusive 
|V..»| measurements which differ by ~ 3.50 (61). |Viv|* is directly proportional to the 
decay rate B° + Xu“ vu, where X, is a meson containing a u quark. Theoretical 
input from Heavy Quark Effective Theory and lattice calculations are also necessary 
to calculate |V,,p|, although several of the theoretical uncertainties cancel in the ratio 
to calculate the side. 


\Viol/|Veo| is a very difficult measurement at LHCb due to presence of a neutrino, 
the identification of which was never in LHCb’s original plans. Although the 
B Factory favoured channel B° > nuv, cannot currently be identified at LHCb, 
the equivalent baryonic channel A, —> p` vy, has been measured. The signal is 
separated from the lower-mass backgrounds, shown in Fig. Using form factors 
from lattice calculations, LHCb measures 


IVa] = (3.27 + 0.15 (exp) + 0.17 (theory) + 0.06 (|Ve|)) x 107° 


This is to be compared to the world average of |V,,5| = (3.94 + 0.36) x 1073 [61]. 
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Figure 3.8: The (pu~) mass distribution in the measurement of |Vub|, showing the contribution 
from A, —> pu~ v, decays and the various backgrounds. 


e The side opposite to a 
The mass difference Am, measured in B? mixing (Fig. 3.3), which is dominated 


by the top-quark loop, provides a measurement of the third side of the triangle, 


Via . . . . . . . Ama 
VA This is proportional to the ratio of mixing frequencies 4/ Am 


are calculated from the lattice with a theoretical error of ~5-10% and systematic 
errors largely cancel in the ratio. 


. Corrections 


Following the measurement by LHCb of the mixing parameters presented above, 
the ratio |Vig/Vis| is 0.210 + 0.001 + 0.008 [49]. Systematic errors can be reduced in 
the future by improved lattice QCD calculations. 
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3.4 Other CP violation measurements 


3.4.1 B, weak mixing phase ¢, in B® > J/wo 


It can be seen in Fig. that V;s appears twice in the B? -p° mixing process, introducing 
a relative “weak mixing phase”, of ¢, to fourth order in À. The B, mixing phase can be 
measured in the channel B? — J/wd, which is governed by a single tree-level diagram 
with a negligible penguin contribution. Hence this mode is the strange-quark analogue of 
the golden mode B® + J/wWK in the B® system. In the B, system CP asymmetry arises 
from the interference of the B? + J/w@ with the mixed process B? > B? — J/wd. In 
the Standard Model, ¢, is expected to be very small, ~ 0.036 + 0.002 rad 52], hence this 
channel is a very sensitive probe for new physics. 

LHCb reconstructs B? — J/Yġ events in the decay modes J/y —> u"u™, and 
o> KK- [75]. This B? final state is an admixture of CP-even and odd contributions, 
therefore an angular analysis of decay products is required. Good tagging performance of 
B? and B? is important, with a total tagging power in this analysis of 4.73 + 0.34. 

The fitted value of ¢, is correlated with AT, the width difference of the B? mass 
eigenstates. The decay B? + J/watm is also added to improve the sensitivity [76]. 
Contours in the (¢,, AI’,) plane are plotted in Fig. The LHCb measurements are 
AT, = 0.0816 + 0.0048 ps~! with the CP-violating phase ¢, = —0.041 + 0.025 rad. 
Figure also shows the world-averaged measurement of ¢, versus AI, showing the 
LHCb result in combination with those from other experiments and the Standard Model 


expectation (61). 
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Figure 3.9: The measurement of the B? weak mixing phase in the (¢,, ATs) plane, showing the 
LHCb result in combination with those from other experiments, the Standard Model expectation, 
and the lo contour of the combination [61]. 
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3.4.2 CP violation in charm 


The Standard Model prediction of CP violation in the charm system is expected to be 
very small ØO (1074) + O (1073), where CP violation can arise in Cabibbo-suppressed 
(CS) decays in the interference between tree and penguin amplitudes. In particular LHCb 
has measured asymmetries in the direct CP-violating channels D°(D’) — mtn and 
D(D°) + K+K- [77]. 

In the LHCb analysis, D° and D decays are identified via two self-tagging decay 
paths. “Prompt” decays (D decays originating from the primary vertex) are characterised 
by the presence of a “soft” low-momentum pion from a D* i.e. D** + Drt and the 
charge-conjugate mode. “Semileptonic” decays are secondary D’s which originate from 
prompt B decays, i.e. BY > D tx and its charge-conjugate state. 

The raw asymmetry (A) for D? + hth~ decays (h = K or r) is defined as 


N(D= )=N(D = f) 
N(D > f)+N(D = f) 
which includes both physics and detector terms: A = Acp + Ap + Ap. Detection 
asymmetry arises from small charge differences associated with the mẹ or “*. Production 
asymmetry arises from different production rates of D* and B in pp collisions. 'To 
eliminate these two contributions and cancel associated systematics, the AAcgp parameter 
is measured in LHCb: 


AAcp = A (KYK?) —A(ata7) = Acp (KTK—) — Acp (ntr). (3.10) 

The raw symmetries are obtained from mass fits, and then by simply counting the numbers 
of D’s decaying to m+” and Kt K—, respectively. 

A measurement performed with Run 1 and Run 2 LHCb data combined gives 


AAcp = (—15.4 + 2.9) x 1074. This is a 5.30 measurement of CP violation in the charm 
system and opens a new window for the study of CP violation in the future. 


A(D > f)= (3.9) 


3.4.3 CP violation in beauty baryons 


CP violation has been observed in B, K, and D decays, but not yet in baryon decays. A 
search for CP violation in the multi-body mode A? + p*a7~atn~ decays was performed 
on LHCb Run 1 data [78]. This decay proceeds via tree and loop diagrams with similar 
contributions and through numerous intermediate resonances, enhancing the possibility 
for CP violation, although in areas where re-scattering effects can play a role. A 3.30 
deviation from CP symmetry was observed, however introducing 6.6 fb~! of Run 2 data 
has not confirmed this result [79]. Hence this measurement awaits further statistics, and 
will be improved when cleaner 2-body B-baryon decays can be added to the study. 


3.4.4 CP violation in charmless 3-body B+ decays 


Yields of Bt > a+ Kt K~ and B© + x K~K* decays show striking asymmetries in the 
region of phase space dominated by re-scattering effects [50]. Similarly huge CP-violating 
effects between Bt > ntntr~ and B7 > n~n ~nt decays are observed in a region of 
phase space including the p°(770) and f2(1270) resonances (31}|82]. Figure [3.10] shows an 
example of the spectacular asymmetries observed in B —> tk K decays, which exceed 
50% at low values of KK- mass. 
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Figure 3.10: Yields of (top) Bt > a+K*K~ and (bottom) B© > 2 K~K* decays as a 
function of K+K mass squared. 
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4 Rare decays 


Within the SM the interplay of weak and Higgs interactions implies that Flavour Changing 
Neutral Currents (FCNCs) can occur only at higher orders in the electroweak interactions 
and are strongly suppressed by the GIM mechanism. This strong suppression makes 
FCNC processes natural candidates to search for physics beyond the SM. If the new 
degrees of freedom do not have the same flavour structure of the quarks/leptons-Higgs 
interactions present in the SM, then they could contribute to FCNCs at a comparable (or 
even larger) level to the SM amplitudes. 

In B-meson decays, experimenters have measured b > s and b > d quark transitions, 
while c + u and s —> d transitions have been measured in D-meson and K-meson decays, 
respectively. At first order, these transitions can occur through two kinds of Feynman 
diagram shown in Fig. |4.1| The first corresponds to the so-called “box” diagram and 
describes the mixing between neutral mesons, discussed in Section [$| the example of 
Fig. [4.1] shows B? mixing. The second kind of diagram, the so-called ” penguin” diagram, 
is responsible for a large variety of FCNC rare decays. The example shown in Fig. |4.l]is 
that of a b— sf*é€ transition. In particular, if the radiated bosons are of electroweak 
type (Z, W or 4-like), the uncertainties in the calculation of the SM predictions due to 
non-perturbative QCD effects are drastically reduced as compared with the case where a 
gluon is radiated. These ”electroweak penguins” are the subject of the discussions in this 
section. 
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Figure 4.1: Examples of loop processes within the SM that allow FCNC b > s quark transition. 
On the left is an example of a box diagram and on the right an example of a penguin diagram. 


Before the first physics run of the LHC accelerator in 2010, the main contributors to the 
study of rare B and D-meson decays were the B Factory experiments (BaBar and Belle) 
and the Tevatron experiments (CDF and D0). However the production rate of bb pairs in 
the ete~ B Factories was typically five orders of magnitude smaller than at the LHC. In 
addition, the lower pp collision energy of the Tevatron (with a correspondingly reduced bb 
cross section which is proportional to the collision energy), and the detector’s reduced 
trigger acceptance for rare B and D-meson decays, implied that LHCb would already be 
the most sensitive experiment even after 1 fb™t accumulated in 2011. For example, prior 
to the LHC, the rarest B-meson decay ever measured was B(Bt > Ktpty-) ~5 x 1077 
with a 20% precision; Belle and BaBar had analyzed O(200) B° > K*°¢té- events and 
CDF had reached a sensitivity of B (B? —> utu) < 3 x 1078 at the 95% C.L. As will be 
clear from the next sections, LHCb has already reached after Run 1 of the LHC an order 
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of magnitude larger statistical power. 


4.1 Bo. > pt pe 


The pure leptonic decays of K, D and B mesons are a particular interesting case of 
electroweak penguins, where the final quark leg in Fig. (right) needs to be swapped to 
the initial state. The helicity configuration of the purely leptonic final state suppresses 
the vector and axial-vector contributions by a factor in terms of masses proportional 


2 
to | -| . Therefore, these decays are particularly sensitive to new (pseudo-) scalar 


interactions. In the case of B° and B°-meson decays, the contribution of the absorptive 
part can be safely neglected. As a consequence, the rate is well predicted theoretically 
: B(B°> pty) = (3.66 + 0.05) x 107° and B(B° > uu) = (1.03 £0.05) x 1071°. In 
the B® case, this prediction corresponds to a flavour-averaged time-integrated measurement, 
taking into account the correction due to the non-negligible width difference. 

The experimental signature is sufficiently clean to reach an expected signal over 
background ratio S/B ~ 3 for the B? decay, assuming the SM branching fraction. The 
main background in the region around the B® invariant mass is due to combinations of 
uncorrelated muons and can be estimated from the mass sidebands. The most important 
handle to reduce this combinatorial background is the invariant mass resolution of the 
experiment, 7 ~ 23 MeV, which is also crucial to differentiate between B° and B® decays 
(Am ~ 87 MeV). Moreover, the large fraction of B—> hth” decays is an important source 
of background due to hadrons being misidentified as muons in the region around the B® 
mass (this background is very small in the B? mass region). 

Given their experimental detector resolution and trigger acceptance during Run 1 
and Run 2, the CMS experiment with 61 fb~’ of data has similar sensitivity to the 
LHCb experiment with 4.4 fb~! of data collected in the same period. Both experiments 
have provided a clear observation of the decay B? —> utu (CMS and LHCb observe 
about 90 and 30 B? signal candidates with a background of about 40 and 10 events 
respectively within the invariant mass interval expected to contain 95% of the signal). 
Neither experiment has yet reached the sensitivity to observe the decay B° > utu”. 
The invariant mass distribution obtained by the LHCb experiment is shown in Fig. 
and the LHCb and CMS results are B(B°—> pt u~) = (3.0 + 0.6493) x 107°? and 
B(B° > ptp-) = (2.9t%7 + 0.2) x 107° respectively, while for the B° decay only limits 
can be quoted. The ATLAS collaboration has also published after a measurement 
using 51 fb~! of data, B(B°—> ptp-) = (2.8778) x 107° in agreement with the previous 
results. 

The B? — pt very rare decay has been searched for ever since the discovery of B 
mesons, around 40 years ago. Thanks to the ingenuity and persistence of the experimenters, 
it has been eventually measured at the LHC and found to be in agreement with the 
SM within current uncertainties, as shown in Fig. Over the next decade it will be 
extremely interesting to see how the measurement of B (B°— pt ~) evolves, for which 
only upper limits are currently available. 
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Figure 4.2: Invariant mass distribution of dimuons selected from Ref. [84]. Superimposed on the 
data points in blue (solid line) is the combined fit, and its components as quoted in the insert. 


4.2 B°> K”? u*u 


The language of effective field theory is used to parameterise NP contributions in terms of 
a sum of local four-fermion operators (Q;) which depend only on SM fermions modulated 
by Wilson coefficients (C;) which in turn depend on the heavy degrees of freedom i.e. NP 
particles. The decay B? — K*° utu" is the so-called “golden mode” to test new vector 
(axial-vector) couplings, i.e. the Cg and Cio Wilson coefficients contributing to the b > s 
transition. The B° > K*° utu channel also complements the b > sy decay which is 
mostly sensitive to NP dipole operators (i.e. Cz). and also the B? — utu decay which is 
mostly sensitive to NP (pseudo-) scalar operators (i.e. Cg and Cp). The charge of the 
pion in the decay K* — Kr defines the flavour of the B meson and an angular analysis 
can be performed unambiguously to test the helicity structure of the electroweak penguin. 

The above system is completely defined by four variables: q?, the square of the invariant 
mass of the dimuon system, 0;, the angle between the positive lepton and the direction 
opposite to the B-meson in the dimuon rest frame, 0g, the equivalent angle of the Kt 
in the K* rest frame and @ the angle between the two planes defined by (K,7) and 
(ut, u~) in the B-meson rest frame. The four-fold differential distribution contains a 
total of eleven angular terms that can be written in terms of seven qg?-dependent complex 
decay amplitudes. These amplitudes can be expressed in terms of five complex Wilson 
coefficients (C's, Cp, C7, Cy and Cio), their five helicity counterparts and six form-factors, 
which play a role of nuisance parameters in the fit. 

The LHCb experiment, with 3 fb~' of data collected in Run 1, has triggered and 
selected about 2400 B? > K*°yu*p- candidates in the range 0.1 < q? < 19 GeV? with 
signal over background (S/B) > 5. This is about one order of magnitude larger than 
the samples available at previous experiments (BaBar, Belle and CDF) and similar to 
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Figure 4.3: Upper limits (B°— utu) and eventual measurements of the B2> ut u~ branching 
ratio by several experiments over the last 40 years. 


the samples collected by ATLAS and CMS with ten times the luminosity, however with 
significantly worse S/B. 

The statistics and the quality of the data accumulated by the LHCb experiment allows 
for a full angular analysis of B? — K*°y*p- decays to be performed for the first time. 
The results of this “tour de force” analysis mostly agree with SM predictions, however 
with some hints of disagreement for some specific distributions. In Fig[4.4] two examples of 
the CP-averaged angular coefficients (i.e. the average of the coefficients measured with B° 
and B® decays) are shown as a function of q?. For these two examples, Arg (modulating 
the sin? 0 x cos 0; angular term) and S; (modulating the sin(20x) x sin 0; x cos ¢ angular 
term) seem to agree less well with SM predictions. However these are early days, and 
more data will be required (the Run 2 data analysis will be released soon). Also a careful 
reassessment of the SM uncertainties are needed before drawing definitive conclusions. 

Several authors have already attempted to see if the overall pattern of the angular 
measurements is consistent with a given value of the relevant Wilson coefficients. As 
previously discussed, the inclusive b + sy measurements strongly constrain non-SM values 
for Cy. The scalar C's and pseudo-scalar C'p coefficients are constrained, for example, by 
the measurement of the branching fraction of the very rare decay B? — y*p~. Therefore, 
the small disagreements observed in the angular analysis of the decay B? > K* u" u 
and other decays, seem to be consistent with a non-SM value of the Cg Wilson coefficient, 
as can be seen in Fig. [4.5] taken from Ref. [88]. 
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Figure 4.4: Two examples of the CP-averaged coefficients in the B? —> K*° u" u7 angular terms 
as a function of q?. The shaded boxes show the SM predictions taken from Ref. [89]. 
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Figure 4.5: Constrains on the contribution of NP to the real parts of Co and Cio at the lo and 
2c level taken from Ref. [88], corresponding to different measurements as indicated in the insert. 
The dotted lines correspond to the status before the latest updates for the 2019 physics winter 
conferences. 


4.3 Lepton Universality 


In the SM, the electroweak couplings of leptons are flavour independent, or lepton 
“universal”. However, this may not necessarily be the case for new particles beyond the 
SM. In particular, if the hints described in the previous section are an indication of new 
particles modifying the penguin diagram in Fig. |4.1| it is interesting to measure the ratios 
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of branching fractions in decays of different lepton families. For example, the ratios 


0 0 
eg > Xun) , we | (Oe mele (44) 


between B decays to final states with muons and electrons, where X is a hadron containing 
an s-quark or d-quark, are predicted to be very close to unity in the SM |90}92]. The 
uncertainties from QED corrections are found to be at the percent level {93}. 

LHCb has measured the above ratio in several channels. Using Run 1 and part of 
Run 2 cane (4.4 fb~'), LHCb measures Rg = 0.846+%:069 (stat) +0- 0i¢(Syst) in the range 
1.1 < q? < 6 GeV’, about 2.50 below the SM prediction, and using only Run 1 data 
(3 fb-') LHCb measures [95] Rg = 0.69*}57(stat) + 0.05(syst) in the same q? range, 
with a similar level of disagreement with the SM prediction. The latest Rx results from 
LHCb in bins of q? can be seen in Fig. compared with previous results from the 
BaBar and Belle collaborations. The consistency between different experiments, and 
different channels, although with very different precision, has motivated many theoretical 
studies that relate these hints of lepton non-universality with the discrepancies described 
in the previous section. Figure shows the status of the compatibility of both sets of 
measurements when assuming that only the bsuu Wilson coefficients are modified (and 
bsee coefficients are as predicted by the SM). Whilst the initial results showed remarkable 
consistency between different sets of measurements as shown by the dotted lines in Fig. |4.5} 
the 2019 latest updates show a less clear picture. As of today, it is difficult to draw reliable 
conclusions, and more data is eagerly awaited. 
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Figure 4.6: Measurement of Rx as a function of the q? bins by different experiments. 
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5 Spectroscopy 


A deep understanding of quantum chromodynamics, the theory of strong interactions, is 
vital for precision tests of the Standard Model and in searches for new physics beyond. 
QCD is intensively tested in deep-inelastic-scattering processes and heavy vector boson 
production, however in the low-energy regime, there is a lack of precise QCD predictions. 
QCD, being a non-pertubative theory, does not calculate hadron properties, namely masses 
and decay widths from first principles. Alternative theoretical approaches are developed, 
such as heavy quark effective theory, heavy quark expansion or lattice calculations. 
These approaches require verification with experiment in various regimes, e.g. testing 
the agreement with data for hadrons with different quark content and quantum numbers. 
Spectroscopic measurements of hadron masses and widths or lifetimes provide a wide 
variety of tests for QCD models. 

Huge production cross sections of charm and beauty in high-energy pp collisions in 
the forward region at LHCb (96}112], together with a good reconstruction efficiency, 
versatile trigger scheme and an excellent momentum and mass resolution, opens up 
exciting opportunities for spectroscopy measurements. The employment of LHCb’s pow- 
erful hadron identification system enables a substantial reduction in the 
combinatorial background specific to high-energy hadron-hadron collisions. The unique 
hadron identification is also especially important for spectroscopy measurements involving 
charged kaons and/or protons in the final state. The excellent momentum and vertex 
resolutions provided by the LHCb tracking system allows unprecedented precision on 
mass and width measurements: indeed the most precise measurements of mass for all 
open beauty particles and lifetimes of all open heavy flavour particles, currently result 
from the LHCb experiment |114|. Additional good control over the momentum scale and 
the detector alignment also allows the natural widths of hadronic resonances to 
be probed with world-leading sub-MeV precision 123}. 

Almost 30 new hadrons have been discovered using the Run 1&2 data-set, among 
them several open charm mesons and baryons, double-charm baryon, a charmonium state, 
several beauty mesons, fifteen beauty baryons, and pentaquark-like P.(4312)*, P.(4400)* 
and P.(4457)* states. The quantum numbers of many hadrons, especially for charm mesons 
and charmonia-like exotic candidates, have been determined using amplitude-analysis 
techniques. In the case of exotic particles such as pentaquark and tetraquark candidates, 
the determination of quantum numbers is vital for the understanding of their nature. 
Large statistics and low levels of background in decays of beauty hadrons with J/Y —> u* u~ 
or W(2S)—> pty final states allow the study of the spectroscopy of light unflavoured or 
strange hadrons (125}{144). Selected spectroscopy highlights are discussed in this chapter. 


5.1 Charm hadrons 


Two complementary methods for the study of spectroscopy of charm hadrons have been 
exploited in LHCb: 


e the study of promptly produced charm hadrons, 


e the study of charm particles produced in weak decays of beauty hadrons, e.g. from 
exclusive B— D“ am decays. 
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The first technique allows the most efficient exploitation of the huge prompt cé production 
cross section in high energy hadron collisions, but this is usually affected by a large 
background from light hadrons produced at the pp collision vertex. The second technique 
exploits a full amplitude analysis of the exclusive decays of beauty hadrons and therefore 
involves much lower statistics, however the method often allows the determination of 
quantum numbers of the charm hadrons. The study of Dz final states enables a search 
for natural spin-parity resonances, (P = (—1)/, labelled as D*) whilst the study of D*r 
final states provides the possibility of studying both natural and unnatural spin-parity 
states, except for the J? = 0+ case, which is forbidden because of angular momentum 
and parity conservation. In inclusive D“)q production, the production of any J? state 
is permitted. An amplitude analysis of B decays allows a full spin-parity analysis of 
the charmed mesons present in the decay. 

Both the above approaches are complementary, and have resulted in discoveries of 
several new charm hadrons, amongst them the excited charm mesons. Many previous 
meson and baryon states, discovered earlier by other experiments, have been confirmed 
with high statistical significance, and their masses and widths have been measured with 
high precision. For many, the quantum numbers were either measured or constrained. 


Excited charm mesons in the D*tz~, Dtr” and D°®x* spectra have been 
studied in LHCb using a 1 fb~! data-set collected at \/s = 7 TeV [145]. Four resonances 
labelled Do(2550), D%(2600), D(2740) and D3(2750) have been observed. The Do(2550), 
and D(2740) decay angular distributions are consistent with an unnatural spin-parity, 
whilst the D*(2600), and D3(2750) states are assigned natural parities. For the Do(2550), 
meson, angular distributions are consistent with a J? = 07 assignment, however for the 
other states, no definite assignment exists. 

The excited charm mesons were also studied in B-decay amplitude analyses of 
B?— D°x+n- and B-—- Dtx-nt decays using the Run 1 data-set and 
Bo > Dtr ort using a data-set corresponding to 4.7fb~', collected at ys = 7, 
8 and 13 TeV. The D3(2760)° and D3(3000)° states were observed for the first time [147], 
and the most precise determination of masses, widths and quantum numbers has been per- 
formed for the Dj(2400)~, D3(2660)~ , D3(2760)~ [146], D3(2660)°, D}(2680)°, D3(2760)°, 
D%(3000)° and for the D,(2420)°, D:(2430)°, Do(2550)°, D*(2600)°, D2(2740)°, 
D3 (2750)° states. 


Excited D} mesons were studied both in inclusive production in the Dt K§ and 
D°K* spectra using a data-set of 1fb~' collected at \/s = 7TeV and through 
the amplitude analysis of B? + D°K~zxt decays using Run 1 data [150]. It has been 
shown that the previously reported D*,(2860)* state is an admixture of 
spin-1 and spin-3 resonances. The masses and width of new states, dubbed D*,(2860)* and 
D*,(2860)*, are precisely measured, as well as the masses and widths of the D*,(2573)* 


and D*,(2700)* [149]]150}. 


Excited Af baryons were studied in their decays to the D°p final state via the am- 
plitude analysis of A? baryon decays using an integrated luminosity sample of 3 fb! 
collected at ys = 7 and 8 TeV (153}. The analysis uses a sample of 11 212 + 126 sig- 
nal A? + D°pr~ decays, where the D° mesons are reconstructed in thek~7z* final state. 
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The amplitude fit is performed in the four-phase space regions in the Dalitz plot. For the 
near-threshold mp», region, an enhancement in the D°p amplitude is studied. The en- 
hancement is consistent with being a resonant state, dubbed the A,(2860)*, with quantum 
numbers J? = ae and with the parity measured relative to that of the A,(2880)* state. 
The other quantum numbers are excluded with a significance greater than 6 standard 
deviations. The phase motion of the component with respect to the non-resonant 
amplitudes is obtained in a model-independent way and is consistent with resonant be- 
havior. The mass of the 4.(2860)* state is consistent with predictions for an orbital 
D-wave At-excitation with quantum numbers F based on the nonrelativistic heavy 
quark-light diquark model and from QCD sum rules in the HQET framework. Also the fit 
allowed the most precise determination of the masses and widths of the known resonances 


A.(2880)* and 4.(2940)*, as well as constraining their quantum numbers. 


Three excited =*° baryons have been observed in the A} K~ mass spectrum using 
the Run 2 data-set [154]. The mass difference 6m = m(AtK7~) — m(At) — m(K7) spec- 
trum for selected At K~ combinations is shown in Fig. |5.1, Three narrow structures, 
denoted £,(2923)°, £.(2939)° and =,(2965)° are clearly visible with a significance ex- 
ceeding 200 for each signal. The data and fit show the least compatibility in the region 
ôm ~ 110 MeV, that could be evidence for a fourth =*° state. Figure [5.1[right) shows 
the dm distribution for the signal samples, where a structure in this region is added into 
the fit. A large improvement in the fit quality is achieved. 
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Figure 5.1: Mass difference Am = m,+;,— — M,+ — mg- distributions 154]. Fits accounting 
for (left) three and (right) four excited =*° states are superimposed. 


Five narrow excited 2* baryons have been observed in the =A mass spectrum 
using the Run 1 data-set (155). A large sample of Æ+ candidates were reconstructed in 
their Cabibbo-suppressed mode =*+— pK~zt. In total around 1.05 x 10° 2+ pk-a«t 
candidates with with a purity of 83% were selected, shown in Fig. [5.2{left). The mass 
distribution for S+ K7 combinations is shown in Fig. [5.2(right), and five narrow peaks are 
clearly visible. The natural widths of the peaks are found to be between 0.8 and 8.7 MeV, 
and two of them, named (2,(3050)° and (2,(3119)°, are found to be extremely narrow, with 
95%CL limits of 1.2 and 2.8 MeV, respectively. It is found that the fit improves if an addi- 
tional broad Breit—Wigner function is included in the 3188 MeV mass region. This broad 
structure may represent a single resonance, be the superposition of several resonances, be 
a feed-down from higher states, or some combination of the above. The interpretation of 
the narrow states is still an open question. The naive quark model expects five states in 
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the region, but some have to be relatively broad. The molecular model predicts two states 


with J? = es and two with J? = and three of the observed states are in remarkable 


agreement, both in mass and width, with this hypothesis (156). 


80000 —— 


400- LHCb -4 
60000} 


40000 - 


Candidates / (1 MeV) 
Candidates / (1 MeV) 
Ww 
© 
S 
| 


20000 L 


[a a a a S E A E ] + 
2440 2460 2480 2500 3000 3100 3200 3300 
m(pK 7+) [MeV] mZtK ) [MeV] 


ao 
— 


Figure 5.2: (left) Distribution of the reconstructed invariant mass Mpg-r+ for all &7 candidates. 
The solid (red) curve shows the result of the fit, and the dashed (blue) line indicates the fitted 
background. (right) Distribution of the reconstructed invariant mass of = K~ combinations; 
the solid (red) curve shows the result of the fit, and the dashed (blue) line indicates the fitted 
background (155). The shaded (red) histogram shows the corresponding mass spectrum from 
the Z+ side-bands and the shaded (light grey) distributions indicate the feed-down from partially 
reconstructed 2,(X)° resonances. 


5.2 Double-charm baryons 


Three weakly decaying states with charm number Č = 2 are expected in the quark model: 
one isospin doublet Z.. and one isospin singlet Rees, each with spin-parity J? = m 
The properties of these baryons have been calculated with a variety of theoretical models. 
In most cases, the masses of the Zee states are predicted to lie in the range 3500 to 
3700 MeV/c? (157). The masses of the =4* and ZF states are expected to differ by only 
a few MeV/c? due to approximate isospin symmetry. Most predictions for the lifetime 
of the =* baryon are in the range 50 to 250 fs, and the lifetime of the =** baryon is 
expected to be three to four times longer at 200 to 700 fs, While both are expected to 
be produced at hadron colliders the longer lifetime of the Ætt baryon should make it 
significantly easier to experimentally observe than the =*, baryon. 

Experimentally, there is a longstanding puzzle in the =,. system. Observations of 
the =* baryon in the At K~z* final state at a mass of 3519 + 2MeV/c? with signal 
yields of 15.9 events over 6.1 + 0.5 events background (6.30 significance), and 5.62 events 
over 1.38 + 0.13 events background in the final state pD* K- (4.80 significance), were 
reported by the SELEX collaboration [158}/159). The SELEX results included a number 
of unexpected features, notably a short lifetime and a large production rate relative to 
that of the singly charmed At baryon. The lifetime was reported to be shorter than 
33 fs at the 90% confidence level, and SELEX concluded that 20% of all At baryons 
observed by the experiment originated from =f, decays, implying a relative =... production 
rate several orders of magnitude larger than theoretical expectations. Searches from 
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the FOCUS [160], BaBar [161], and Belle experiments did not find evidence for 
a state with the properties reported by the SELEX collaboration, and neither did a search 
at LHCb with data corresponding to an integrated luminosity of 0.65 fb~! (163). However, 
because the production environments at all the above experiments differ from that of 
SELEX, which studied collisions of a hyperon beam on fixed nuclear targets, these null 
results do not exclude the original observations. 

LHCb has searched for the Z4" decaying into At K~atxt using a sample of pp col- 
lision data at 13 TeV, corresponding to an integrated luminosity of 1.7fb~’. A highly 
significant structure is observed in the mass spectrum, where the At baryon is recon- 
structed into pKr", shown in Fig. [5.3{ left). The structure is consistent with originating 
from a weakly decaying particle, identified as the doubly charmed baryon =7*. The 
observation of the state is confirmed using an additional sample of data collected at 8 TeV. 
Soon after, the observation was further confirmed by observing the same state in the decay 
Sit + Sir, shown in Fig. [5.3{(right). The mass of the 5+ state was measured to 


cc 
be in very good agreement with the value measured in the =t*+ + AtK-atnt decay 
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channel |164|. The lifetime and mass of the 5** baryon were also precisely mea- 
165/166 


sured , where the lifetime favours smaller values in the range of the theoretical 
predictions. 
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Figure 5.3: (Left) Mass distribution of At K~a*x* candidates with fit projections overlaid : 
(Right) Mass distribution of S/t candidates [164]. 


5.3 Beauty hadrons 


Excited B* and B® mesons have been investigated in the mass distributions of Btr- 
and B°r+ combinations using a 3fb~’ data sample at 7 and 8TeV. The Bt and B® can- 
didates were reconstructed through the Bt — D°xt+, Bt = D°xtata-, Bt = J/WKt, 
B? + Dont, B} > Dontrtr and B° > J/Wk* decay chains. Samples of about 
1.2 million B° and 2.5 million B+ candidates have been obtained with purity depending 
on decay mode, but always better than 80%. The Btz~ and B°r* mass spectra with 
requirements that pr > 2 GeV/c are shown in Fig. |5.4| where ten peaking structures are 
reconstructed. Out of these, six narrow low-mass structures correspond to the decays 
of the four B,(5721)°+ and B3(5747)°* states observed by the CDF and D0 collabora- 
tions Tal B,(5721) > B*r and B2(5747) + B®r. The high statistics of LHCb 
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has allowed the most precise measurements of the masses and widths of the B,(5721)°* 
and B}(5747)°* states to be made. 

In addition to the six low-mass structures, four wider high-mass structures are observed, 
particularly prominent at high pion transverse momentum. These structures are consistent 
with the presence of four new excited B mesons, labeled B;(5840)°* and B,(5960)°*, 
whose masses and widths are obtained under different hypotheses of their quantum 


numbers (171). 
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Figure 5.4: Spectra of Am = Mgr — mp —m, for (left) Btr and (right) B°r+ candidates 


with pr(7) > 2 GeV/c [171]. 


Orbitally excited B? mesons have been studied using only 1 fb~! of data, collected 
at \/s = 7 TeV. The Bt K~ mass spectra were investigated with the Bt mesons being 
reconstructed in four decay modes. Previously, two narrow peaks had been observed in 
the BK- mass distribution by the CDF and DO collaborations [172/173], named the 
B* (5830)° and B*,(5840)°. They are putatively identified as members of jg = 2 HQET dou- 
blet . The two states are also visible in LHCb data, here as three narrow peaks shown 
in eset, corresponding to the decays B,;(5830)° + B*t K~, B*,(5840)° > B*t K7, 
and B*,(5840)° > Bt K—, where a soft photon from B*t + Bty is undetected. This is 
the first observation of the B*,(5840)° + B*t K- decay mode and a J? = 2+ assignment 
is favoured for this state. Large statistics, low background and LHCb’s excellent mass 
resolution has allowed the first determination of the B*,(5840)° width as well as the most 
precise mass measurements of both states. Due to the small energy release, the position 
and the shape of the peaks depends on the mass of the B** states, allowing the most 
precise determination of the mp.+ mass, as well as the mass difference mg:+ — mp+. 


Excited B+ mesons have been searched for via their decays into the Bia*a~ final 
state. A wide peak, interpreted as the B.(29)*, was observed by the ATLAS collaboration 
using a sample of about 300 reconstructed Bt => J/wat candidates [175], with a large 
relative production rate with respect to the base B} state. LHCb has searched for this 
state using 2fb~' of data collected at ys = 8 TeV, observing a sample of 3325 + 73 
reconstructed B7 — J/wat decays. No signal is observed and an upper limit on the 
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Figure 5.5: (left) The mass difference distribution mpg+g- — mg+ — mg- in the 


B+ K~ mass spectra. The three peaks are identified as (left-to-right) B,1(5830)° > B*t K7, 
B*,(5840)° + B**K~, and B*,(5840)° + Bt K~—. The total fit function is shown as a solid blue 
line, while the shaded red region is the spectrum of like-charge Bt Kt combinations. The inset 
shows an expanded view of the Bs1(5830)°/B*,(5840)° + B**K~ region. (right) Distribution of 
AM = Mpt+7+_- — Mpy in the Bian mass spectra with fit results overlaid. The same-sign 
distribution has been normalised to the data in the B, (2.5) side band region. 


relative production rate has been obtained [176]. This upper limit is smaller than the 
relative production rate reported by ATLAS. 

In 2018, the CMS collaboration, using a huge data set corresponding to 143 fb~' 
collected at \/s = 13 TeV and containing of 7629+ 225 signal Bt > J/wx* decays, reported 
observation of a doublet of two narrow states, interpreted as spin-triplet B.(2S)** and spin- 
singlet B.(2S)* states (177]. These observations were confirmed by LHCb using a 8.5 fb”! 
data-set collected at \/s = 7,8 and 13 TeV [178], with 3785 +73 signal B} + J/yr* decays. 
Two narrow peaks with a width compatible with the detector resolution are seen in 
the mass-difference M p++,- — Mpy spectrum, shown in Fig. ee The local (global) 
significances of the two peaks are estimated to be 6.80 (6.30) and 3.20 (2.20) for the 
low-mass and high-mass states, respectively. The low-mass signal is interpreted as the 
spin-triplet state Be(2S)*t, decaying into Bžtr*mn with the subsequent decay of the B** 
into By. The high-mass peak is attributed to the decay of the spin-singlet B,(2S)* state 
into the Brr” final state. 


Excited A? baryons were discovered in the A?r+7~ mass spectrum using only 1 fb~! of 
LHCb data, accumulated at ys = 7 TeV in 2011 |117|. The A? candidates were recon- 
structed via A? Ata followed by At — pKr". In total (70.54 + 0.33) x 10° signal 
Ap decays were selected and then combined with rtr pairs. The A?a*a~ spectrum in 
the region 5.90 < MAntat < 9.95 GeV/c? is shown in Fig. |5.6(left), where two narrow 
peaks, consistent with detector resolution, are visible. The significances of the observations 
are 5.2 and 10.2 standard deviations for the low-mass and high-mass peaks, respectively. 
The observed states are interpreted as the doublet of orbitally-excited A,(1P)° states 
with quantum numbers J? = ~ nd 3™. In 2020 the analysis was updated using the full 
Run 1 and 2 data-sets, shown in Fig. right). The masses of these states are measured 
with unprecedented precision (121). 
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Figure 5.6: (Left) Mass spectrum of A?a*m~ combinations (117]. The points with error bars are 
the data, the solid line is the result of a fit, and the dashed line is the background contribution. 
(Right) Mass spectrum of (top) APatm—, (middle) APrtat and (bottom) Afr- a7 combinations 
in the full Run 1 and 2 data-sets [121]. 


Using the full Run 1 and 2 data-sets, the mass spectrum of A?a*+z* combinations, where 
Ap— AFT and AP J/bpK—, was explored for higher masses of A?a*z~ combinations. 
A significant broad structure is found at m ~ 6.150 GeV/c? with a width around 10 MeV, 
shown in Fig. [5. 7[ left) [179]. The mass and width agree well when measured in the 
two decay modes A? Afm and A9 —> J/pK—, where the significance exceeds 26 and 
9 standards deviations, respectively. Since the mass of the new structure is above 
the StF kinematic thresholds, the Artnr mass spectrum is investigated in the 


Art mass regions populated by the oF ag resonances. 

The data are split into three non overlapping regions: candidates with a A?7* mass 
within the natural width of the known XF mass, candidates with a A?a* mass within 
the natural width of the known ¥;* mass, and the remaining nonresonant (NR) region. 
The A?a*z~ mass spectra in these three regions are shown in Fig. [5.7{right). The spectra 
in the X, and X% regions look different and suggest the presence of two narrow peaks with 
very similar widths. The two-signal hypothesis is favoured over the single-signal hypothesis 
with a statistical significance exceeding seven standard deviations. The masses of the two 
states measured are consistent with predictions for the doublet of A,(1D)° states with 
quantum numbers J? = and a 

In 2020, the fifth excited A, state was observed in the APa*n~ mass spectra, using 
the full LHCb Run 1 and 2 data-sets. Two decay modes of the A? baryon were used, 
Ap Atn~ and AP> J/pK~—, and the significance of the new state, denoted A;*”, is 
in excess of 14 and 7 standard deviations in the two decay modes, respectively. This 
is shown in Fig. Unlike the previously observed four narrow A, states, A,(5912)°, 
A,(5920)°, A,(6146)° and A,(6152)°, the new state is rather broad, I = 72 +11 +2 MeV. 
The measured mass and width agree with the interpretation of this state being the first 
radial excitation, the A,(2S)° resonance |180|. This resonance is also consistent with 
a broad excess of events in the A?a*+z~ mass spectrum, previously reported by the CMS 
collaboration [181]. 
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Figure 5.7: (left) The mass distribution of selected A?a+m~ candidates for (top) the A?— Ata 
and (bottom) the A? > J/ypK~ decay modes. (right) Mass distributions of selected Artnr can- 
didates for the three regions in A?z* mass: (top) Xp, (middle) X% and (bottom) the nonreso- 
nant (NR) region [179]. 
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Figure 5.8: Mass distribution of selected (top) A?a*a~, (middle) Artr* and (bottom) Alr- a 
candidates for the (left) A? Afa~ and (right) A? > J/ppK- decay modes 121|. 


Excited s- baryons have been studied in the A9r* mass spectra using the Run 1 
LHCb data-set [182]. In total (234.27 +0.90) x 10° signal A? baryons were reconstructed in 
the decay mode A? — Atz~. Distributions of the energy release in the decay, Q = m Aone = 
MA9 — Mp, are shown in Fig. 5.9} At low values of Q, there are previously-known signals 


from Seg states, observed and characterised by the CDF collaboration 18341184). New 


peaks in the Alr (Art) spectra are visible at Q = 338.8 + 1.7 MeV (336.6 + 1.7 MeV), 
with a local significance of 12.70 (12.60), based on the differences in log-likelihoods between 
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the fits with zero signal and the nominal fit. In the heavy-quark limit, five 1,(1P) states 
are expected, and several predictions of their masses have been made. Since the expected 
density of baryon states is high, it cannot be excluded that the new observed structures 
are the superposition of more than one (near-)degenerate state. Taking into account that 
the predicted mass and width depend on the as-yet-unknown spin and parity, the newly 
observed structures are compatible with being ¥,(1P)* excitations. Other interpretations, 
such as molecular states, may also be possible. 
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Figure 5.9: The distributions of energy release Q = m Azt — MA9 — Mx for selected A candi- 


dates |182|. The points show experimental data. The left (right) plot shows A}m* (A977) combi- 
nations. 


—/— 


=, and =, paryons have been observed in the Sfr mass spectrum using the 
LHCb Run 1 data-set . Signal =? candidates were reconstructed in the final state 
Etr with =. pK i “te peaks are clearly visible in the ôm = Mazor- — Mz — Mr- 
spectrum, shown in Fig. [5.10{left), a narrow state at ôm ~ 3.7 MeV/c? and a broade state 
at ôm ~ 24 MeV/c. No structure is observed in the wrong-sign sample, nor in the =? mass 
side bands. The fitted natural width of the lower-mass state is found to be consistent with 
zero. The fitted yields of the lower and higher-mass peaks are 121 +12 and 237+ 24 events, 
respectively, with statistical significance in excess of 10 standard deviations. The non- 
zero value of the natural width of the higher-mass state, [ = 1.65 + 0.31 = 0.10 MeV is 
significantly different from zero. The signals are interpreted as £; and =}~ baryons. 


The Æ% baryon was first observed at CMS Te later studied in detail by 
the LHCb collaboration using the Run 1 data-set |119|. The 5%° candidates have 
been reconstructed in the decay =}? 5, at, with 5, 3 S0n* ence Z> pK Krt. 
The ôm distribution, defined as Man qt — Ms- — Mrt, is shown} in Fig. |5.10{right). A nar- 
row peak is clearly visible with a fitted signal yield of 232 + 19 events. The non-zero 
value of the natural width of the peak, [ = 0.90 + 0.16 MeV, is also highly significant; 
the change in log-likelihood when the width is fixed to zero exceeds 30 units. No other 
statistically-significant structures are seen. The peak position and the width are consistent 
with, and about a factor of ten more precise than, the CMS measurements [185]. The mea- 
sured width of the state is in line with theory expectations; a calculation based on lattice 
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Figure 5.10: (left) The distribution of the mass difference, dm = Mzor- — Mz0— Ma, 
for ca candidates [118]. The points with error bars show right-sign candidates in the =, 
mass signal region, and the hatched histogram shows wrong-sign candidates with the same 
selection. The curve shows the nominal fit to the right-sign candidates. Inset: detail of the region 
2.0-5.5 MeV/c?. (right) Distribution of dm and the fit for 5; r candidates [119]. 


QCD predicts a width of 0.51 + 0.16 MeV [186], and another using the 3P) model obtains 
a value of 0.85 MeV [187]. The measured production ratio with respect to the =, state is 
measured to be (28+3+1)%, suggesting that in high energy pp collisions at 7 and 8 TeV, 
a large fraction of =, baryons are produced through feed-down from higher-mass states. 


A high-mass excited =, baryon has been observed in the A?A~ and =?n~ mass 
spectra, using a 3.5 fb’ LHCb data-set at \/s = 7, 8 and 13 TeV (188). The A? baryons were 
reconstructed via A? — Ar~ and A9 — Atu-X, with the =? decaying to Z}? —> Stu X, 
with At, £.— pK~7x*. Full reconstruction of the AP baryon allows the excellent resolution 
of the A? K- mass spectra to be reached. In addition, partial reconstruction of A? and 
=? baryons in their semileptonic modes allows a significant increase in the sample of A? 
and =? baryons, where the missing neutrinos do not prevent a peaking structure in the 
spectra of mass differences (m» jo» g — M»,9») and (M» sory- — M»so). The resolution is 
improved by applying the 4-vector constraint (P+ + Pp- + Pmiss)? = Mio, where H7 
stands for At and Æ+ and H? stands for AP and =?. 

The mass-difference spectra are shown in Fig. where the peak locations for 
all three modes are seen to agree well. The statistical significances of the new excited 
baryon, dubbed the =,(6227)~, are found to be 7.90 for the =,(6227)~ > A} K- followed 
by AP Atr, 250 for the 5,(6227)- + A? K- followed by A? —> Atu- X and 7.20 for 
the 5,(6227)~ + =?a- followed by 2? —> Słu X. 


Four narrow excited p states have been observed in the =?K~ mass spectrum 
using the full Run 1 and 2 LHCb data-sets [189]. The Æ? candidates were reconstructed 
in Æt" final states with =.— pKr". After multivariate selection, a low-background 
sample of (19.2 + 0.2) x 103 S?— S27 decays has been obtained. The mass-difference 
spectrum m(=?K~) — m( Z0) for ZLK- combinations exhibits four narrow peaks, shown 
in Fig. The natural widths of the three lower mass states are consistent with zero, 
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Figure 5.11: Spectra of mass differences for 2% candidates, reconstructed in the final states 
(left) APA, with AP + Ata, (middle) APK~, with A? —> Atu X, and (right) =?Pa7, with 
Z0 > Z}u X, along with fits to the data 188}. The top row is for 7 and 8 TeV data 
and the bottom is for 13 TeV. The symbol M* represents the mass after the 4-vector con- 


straint (Put + Pp- + Pmiss)? = Mino is applied. 


while the width of the high-mass state is found to be 1.4*j3 + 0.1 MeV. The peaks have 
local significances that range from 3.6 to 7.2 standard deviations. After accounting for 
the look-elsewhere effect, the significances of the two low-mass peaks are reduced to 2.10 
and 2.60, respectively, whilst the two higher-mass peaks exceed 5c. The observed =? K- 
peaks seen here are similar to those observed in the £} K7 invariant mass spectrum |155). 
Arguably, the simplest interpretation is that the peaks correspond to excited (2, states, 
in particular the L = 1 angular momentum excitation of the ground state, or possibly an 
n = 2 radial excitation. 

Many of the quark-model calculations predict L = 1 states in this mass region, and at 
least some of the states should be narrow. In particular, using the 3P) model, five states 
in this mass region are predicted, with approximately 8 MeV mass splittings; the four 
lightest have partial widths, [(=?K—), below 1 MeV, whilst the one with the largest 
mass has ['(=?K~) = 1.49MeV. Conversely, predictions using the chiral quark model 
indicate that the JP = (3/2)~ and (5/2)~ states are narrow, but the (1/2)~ states 
are wide {190}. Quark-diquark models have also predicted several excited (2; states in 
the region around 6.3 GeV, with mass splittings similar to that observed here, however, 
there are no predictions for the decay widths. Molecular models have also been employed, 
where two narrow J? = (1/2)~ states are predicted at 6405 MeV and 6465 MeV (156), 
however do not match well with the LHCb measurements. 

An alternate interpretation for one or more of the observed peaks is that they 
arise from the decay of a higher-mass excited Q;*~ state to Z (—> =?7°)K—, where 
the 7° meson is undetected. If the mass of a not-yet-observed =/° state is in the region 
Mz + Mro < Mszp < Mz, each of the observed narrow peaks can be interpreted as having 
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come from the above decay, provided that the corresponding excited {2;*~ state is narrow, 
lo- < 1MeV. In this case, their masses can be evaluated as Mo- = Mszp + dMpeak, 
where OMpeak is a measured position of the peak in the m(=?K—) — m(P) spectrum. 
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Figure 5.12: The mass difference 6m = m(=?PK~) — m( 29) for selected =? K~ candidates [189]. 


5.4 Conventional charmonia and bottomonia 


Charmonium states in DD mass spectra near threshold have been studied using 
the full LHCb statistical sample collected in Runs 1 and 2 (191). D? and D* candidates 
were reconstructed in the D? > K-zt and Dt + K-rtr* decay modes. In total, 3.6 x 10° 
D°D®° and 2.0 x 10° D*+D~ pairs have been selected. The mass spectra for D°D° and 
D*D~ combinations are shown in Fig. left), with the zoomed-in region between 
3.80 < Mpp < 3.88 GeV/c? presented in Fig. right). 

Four peaking structures are observed in the spectra. Two of the peaks correspond to the 
known w(3770) and x-2(3930) charmonium states. A narrow peak close to threshold rep- 
resents partially reconstructed xe (3872) > D*°D® decays, subsequently with D*°— D°y 
or D*°— D°r° with the y or 7° meson missing. The narrow peak with mass around 
3840 MeV/c? is identified as a new charmonium state. Its mass value and small natural 
width suggest an interpretation as the 73(1°D3) charmonium state with quantum numbers 
ae a (192). In addition, prompt hadroproduction of the y.2(3930) and w(3770) char- 
monium states have been observed for the first time, and precise measurements of their 
resonance parameters have been performed. 


Observation of the decays Xe — J/pptu- using a 4.9fb~* data-set collected at 
Vs = 7, 8 and 13 TeV enabled the most precise direct determination of the masses of 
the Xa and x,2 states and the width of the Xe2 to be performed with unprecedented 
precision (120). The observation of these decay modes provides opportunity for the precise 
measurements of the x-1,-2 production and polarization, that in turn is vital for tests of 
QCD models of charmonia production. 
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Figure 5.13: (left) The DD mass spectra. The open (red) histogram shows D°D° combinations 
whilst the dashed (blue) histogram shows D+ D~. (right) The mass spectra of (top) D°D®° and 
(bottom) D+ D~ combinations in the narrow mass region 3.80 < mpp < 3.88 GeV/c’, with fits 
superimposed. Different components employed in the fit are indicated in the legend. 


Observation of w2(3823) > J/pata- in Bt (y2(3823) > Jat) Kt decays 
using the full LHCb Run 1 and 2 data-sets, has allowed the most precise determination 
of the mass of the tensor w2(3823) state and the best constrained upper limit of its 
width [123]. The observed mass distribution is shown in Fig. Within the factor- 
ization approach, the branching fraction for the decay B* — y2(3823)K* vanishes, and 
a non-zero value for this branching fraction allows an evaluation of the contribution of 
the Dot poe rescattering amplitudes in the Bt > ccK* decays. 
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Figure 5.14: Distributions of the (left) Jibata K" and (right) J/pntr mass for selected 
Bt — (2(3823) > J/fpata—) Kt candidates |123]. 


5.5 Pentaquarks 


The LHCb collaboration studied A? > J/ppK decays using the Run 1 data-set |193]. 
In total (26.0+0.2) x 10° signal A? candidates were selected, and an anomalous peak in the 
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Japp mass spectrum is observed, shown in Fig. [5.15[ left). If the peak structure represents 
a resonance which strongly decays into J/wp, the minimal valence quarks would be céwub, 
a charmonium pentaquark state. A full six-dimensional amplitude fit with resonance 
invariant masses, three helicity angles and two differences between decay planes has been 
applied to describe the data. The amplitude model in the fit contains 14 well-defined 
A* states and two pentaquark states, labelled as P.(4380)* and P.(4450)+. The projections 
of the fit are shown in Fig. |5.15| The masses and widths of the wide P.(4380)* and 
narrow P.(4450)* states have been measured. The preferred spin-parity assignments are 


(3 : ; J (3 | ; 2 ) and Ge a5), where the first value is the J? assignment given by 
the best fit. 


—s— data 
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Figure 5.15: The (left) mj, and (right) mgp mass distributions, showing superimposed fit 
projections as solid (red) points. The solid (red) histogram shows the background distribution, 
the blue shaded histogram represents the P.(4450)* state, and the shaded purple histogram 
represents the P.(4380)*. Each A* component is also shown. 


Following the first observation, the exotic hadronic character of the J/wp structure 
around 4450 MeV/c? was confirmed in a model-independent way |141|. This analysis gave 
similar results and excluded that the data could be described by the pK” constributions 
alone. Further confirmation comes from the amplitude analysis of the Cabibbo-suppressed 
decay AP J/bpx- [194], where 1885 + 50 signal candidates were investigated. There are 
different theoretical interpretations suggested, including a tightly bound duucé state 
a loosely bound molecular baryon-meson state or a triangle-diagram processes. 

A partial update of the above analysis was made using the full Run 1&2 data sam- 
ple (195). A nine-fold increase of statistics is achieved due to the larger data sample, 
an improved selection criteria, and increased pp— bb cross-section at ys = 13 TeV in 
Run 2. For candidates with a mass consistent with the nominal A? baryon mass, the J/wp 
and pk” mass spectra were investigated. In the distribution of J/p mass, the previously 
reported peaking structure around 4450 MeV/c? was confirmed, and a new narrow peak 
with mass around 4312 MeV/c? was found. The 4*— pK- contributions are clearly seen 
in the Dalitz plot, shown in Fig. (left). 

Since the newly observed peaks are narrow, the full amplitude analysis faces computa- 
tional challenges. This is because resolution effects should be included in the formalism 
which complicates the fitting procedure. Conversely, narrow peaks can not be due to 
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Figure 5.16: (left) Dalitz plot of A? + J/ypK~ candidates. The vertical bands correspond to A* 
resonances. The horizontal bands correspond to the P,(4312)+, P.(4440)*, and P.(4457)* struc- 
tures at Manyp = 18.6, 19.7, and 19.9 GeV’, respectively. (right) Projection of the A*-suppressed 
M Jyp distribution showing a fit with three Breit—Wigner amplitudes and a sixth-order polyno- 
mial background. The mass thresholds for the X.D? and X7.D*° final states are superimposed. 


reflections from A* states, motivating the validity of the one-dimensional fit approach to 
study the J/wp invariant mass. The J/wp mass in the narrow-resonance region together 
with the result of the fit is shown in Fig. [5.16] (right). The previously reported peak 
around the 4450 MeV/c? mass is now resolved into a two-peak structure of P.(4440)+ and 
P.(4457)* states. In total, three narrow pentaquark states are observed. The statistical 
significance of the two-peak interpretation of the previously-reported single P.(4450)* 
structure is 5.40. The statistical significance of a new P.(4312)* state is 7.30. The masses 
and widths of the pentaquark candidates are measured. Taking into account systematic 
uncertainties, the widths are consistent with the mass resolution. Hence, upper limits on 
the natural widths at the 95% confidence level (CL) are obtained. 

In summary, whilst the existence of pentaquark-like resonances is certainly beyond 
doubt, their exact nature is still unclear. They can be genuine five-quark bound states, or 
e.g. near-threshold meson-baryon molecules. More studies are required to clarify this. 


5.6 Charmonium-like exotic states 


The enigmatic X (3872) particle was discovered in B* decays by the Belle bollab- 
oration [196]. Subsequently, its existence has been confirmed by several other experi- 
ments aTi The nature of this state is rather unclear. Among the open possibilities 
are conventional charmonium and exotic states such as D*°D° molecules (200), tetra- 
quarks (201), or their mixtures (202). Determination of the JP? quantum numbers is 
important to shed light on this ambiguity. The C-parity of the state is positive since 
the X (3872) > J/wy decay has been observed [203] 204). The CDF experiment anal- 
ysed three-dimensional angular correlations in a relatively high-background sample of 
2292 + 113 inclusively-reconstructed X (3872) > J/ata-, J/p —> utp decays, domi- 
nated by prompt production in pp collisions. The unknown polarisation of the X (3872) 
limits the sensitivity of the measurement of JPO [205]. A x? fit of JPC hypotheses to the 
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binned three-dimensional distribution of the J/) and n*m helicity angles and the angle 
between their decay planes |206}/208], excluded all spin-parity assignments except for 1** 
and 277. 

Using ys = 7 TeV pp collision data corresponding to 1 fb~' collected in 2011, the LHCb 
collaboration performed the first analysis of the complete five-dimensional angular corre- 
lations of the Bt + X (3872)K+t, X (3872) > J/bata-, Jip > tp decay chain [209]. 
About 38000 candidates passed the multivariate selection in a +20 range around the Bt 
M Jhont+n-K- Mass distribution, with a signal purity of 89%. The Am = m gryata- — Msp 
distribution is shown in Fig. |5.17[left). The fit yields 5642 + 76 and 313 + 26 candidates 
for (2S) > J/batn- and X(3872) > J/prtr signals, respectively. 

The angular correlations in the B* decay carry information about the X (3872) quantum 
numbers. To discriminate between the 1**+ and 2~* assignments, a likelihood-ratio test 
is used. A test statistic t is defined as —2In[£(2~*)/L(1**)]. Positive (negative) values 
of the test statistic for the data, taata, favour the 1**+ (27+) hypothesis. The value of 
the test statistic observed in the data is taata = +99, thus favouring the 1** hypothesis. 
A rejection of the 2~* hypothesis with greater than 5ø significance is demonstrated using 
a large sample of pseudo-experiments. As shown in Fig. [5.17(right ) the distribution 
of t is reasonably well approximated by a Gaussian function. Based on the mean and 
r.m.s. spread of the t distribution for the 2~+ experiments, this hypothesis is rejected 
with a significance of 8.40. Hence the obtained results correspond to an unambiguous 
assignment of the X (3872) state to be 177. 

The above result rules out the explanation of the X (3872) as a conventional 
ne2(1! D2) state. Among the remaining possibilities are the xa (2°P,) charmonium state, 
and unconventional explanations such as a D*°D° molecule tetraquark state charmoni- 
um-molecule mixture. 
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Figure 5.17: (left) Distribution of Am = m(J/hpata) — m(J/p) for Bt > J/pn nK candi- 
dates. The fits of the 7(2S) and X (3872) signals are displayed. The solid blue, dashed red, and 
dotted green lines represent the total fit, signal component, and background component, respec- 


tively. (right) Distribution of the test statistic t = —2 ln ai for the simulated experiments 


with JPO = 27+ (black circles on the left) and with JPC = 1+* (red triangles on the right). 
The value of the test statistic for the data, tgata, is shown by the solid vertical line. 


With a larger 3fb~* data-set at \/s = 7 and 8 TeV, the analysis has been repeated in 
the decay X (3872) > J/pp? without an assumption on the orbital angular momentum 


54 


210|. The analysis confirmed the J?° = 1+* assignment for the X (3872) state and also 
set an upper limit of 4% at 90% C.L. on the D-wave contribution. 

A precise determination of the mass and width of the X (3872) state was performed 
using two minimally overlapping data-sets. The first was the 3fb-' Run 1 data-set 
in which the X (3872) particles were now selected from decays of hadrons containing 
b quarks [122]. The second was the full Run 1 and 2 data-sets using a sample of 
(547.8 + 0.8) x 103 Bt + Jhpntn” K* decays (123}. In both cases the X (3872) was 
reconstructed in the X (3872) > J/wta final state. The mass and width were determined 
from a fit to the J/yx*n~ mass distribution assuming a Breit-Wigner line shape for 
the X (3872) state, measured as follows: 


myx(3872) = 3871.70 + 0.07 +0.07 MeV/c?, Taw = 1.3940.2440.10 MeV, 
Mx (3872) = 3871.60 + 0.06 +0.03 MeV/c?, Tpw = 0.961933 + 0.21 MeV, 


where the first and the second lines correspond to Refs. and [123], respectively. 
The above measurements represent the most precise determination of the mass of 
the X (3872) state and the first measurements of its width. The measured mass corresponds 
to the binding energy ôE, defined as mpoc? + mpx0c? — mx (3872)c”, which is 70 + 120 keV. 

Whilst the proximity of the measured mass of the X (3872) to the D°D*° thresh- 
old favours the interpretation of this state as a D°D*° molecule, the large produc- 
tion cross-section of the X (3872) disfavours this. The pure molecu- 
lar interpretation is further disfavoured by the observation of the decay X (3872) — y (25)y. 
The ratio of the decay rates to y(2S)y and J/w7 final states is very sensitive to the nature 
of the X (3872) state. This is predicted to be in the range (3—4) x 107? for a DD* molecule 
1.2 — 15 for a pure charmonium state and 0.5 — 5 for a molecule-charmonium mixture. 

LHCb performed the search for X (3872) > y(2S)y decays using the decay chain 
Bt - X(3872)Kt+t,X (3872) > y(2S)y. The analysis was based on a 1 fb~* data sample 
at 7TeV and 2fb~* at 8TeV. The significance of the B+— (X (3872) (2S)y) Kt 
signal is 4.4 standard deviations. The branching fraction, normalised to that of the 
X (3842) > J/wy decay mode, is measured to be 


B(X (3842) > o(2S)y) 


= 2.46 + 0.64 + 0.29. 
B(X (3842) > Jhb) 


This result is compatible with, but more precise than, previous measurements |204\|218}, 
and strongly disfavours a pure molecular interpretation of the X (3872) state. 


Structures in the J/~d system have acquired great experimental and theoretical 
interest since the CDF collaboration reported 3.80 evidence (14 + 5 events) for a nar- 
row ( = 11.7183 + 3.7 MeV) near-threshold X (4140) mass peak in a sample of 75+ 11 re- 
constructed Bt + J/poK* decays [219]. Much larger widths are expected for charmonia 
states at this mass, therefore its possible interpretations as a molecular state, a tetraquark 
state, a hybrid state or a rescattering effect have been discussed. The X (4140) structure 
was confirmed by CMS and DO [221222], however searches in Bt > J/poK* decays 
were negative in the Belle 223)[224] and BaBar experiments. 

Using a 0.37 fb! data-set at \/s = 7 TeV (346+20 signal Bt + J/poK* decays) LHCb 
initially found no evidence for the narrow X (4140) structure [226], in 2.40 disagreement 
with the measurement by CDF, seen in Fig [5.18{left). However, using a significantly 
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larger sample of 4286 + 151 Bt — J/~@Kt decays (the Run 1 data-set), with roughly 
uniform efficiency across the entire J/W@ mass region, LHCb performed a full amplitude 
analysis, including resonant contributions from K* resonances decaying into ¢A* and 
possible resonances in the J/Yo system (142]/143). Four resonance contributions labeled 
as X (4140), X (4274), X (4500) and X (4700) with quantum numbers 1**, 1*+, 0** and 
Ott, respectively, are observed, shown in Fig. [5.18{right). The statistical significance 
varies from 5.6 to 8.40. The widths of the states are found to be between 56 and 120 MeV, 
significantly exceeding the narrow-width of the X (4140) reported by CDF. 
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Figure 5.18: (left) The distribution of mass difference m(J/bd) — m(J/p) in a 0.37 fb! data-set 
for selected Bt > J/bdKt candidates [226]. A fit of the X (4140) signal on top of a smooth 
background is superimposed (solid red line). The dashed blue (dotted blue) lines illustrates 
the expected X (4140) (X(4274)) signal yield from the CDF measurement [219]. The top and 
bottom plots differ by the background function used in the fit: (top) an efficiency-corrected 
three-body phase-space; (bottom) a quadratic function multiplied by the efficiency-corrected 
three-body phase-space factor. (right) The distribution of J/po mass in the full Run 1 data-set 
for Bt > J/pġK* candidates (black data points) compared with the results of the amplitude 
fit containing eight K** — Kt contributions (shown with open symbols) and five X > J/yọ 
resonance contributions (shown as hatched histograms) (142]/143}. The total fit is given by 
the red histogram. 


The charged charmonium-like state Z.(4430)~ was first observed by the Belle col- 
laboration in the ~(2S)m~ mass spectrum of B° > Y(2S) K-n decays [227]. The state ap- 
peared as as a narrow (T = 44+17+30 MeV) structure with a significance of 6.50. Later, the 
collaboration performed a full amplitude analysis of 2010 + 50 + 40 B? > Y(2S) K-T" sig- 
nal decays, determining the quantum numbers as J? = 1+, and finding a much broader 
width of T = 2007%5°38 MeV [228]. 

The LHCb experiment has collected about 25, 000 signal B° decays (the Run 1 data set) 
and observed the Z,(4430)~ with a significance exceeding 13.90 (137]. Model-independent 
as well as full amplitude analyses were performed. The spin-parity is confirmed as 17, 
other hypotheses are excluded by at least 9.70. 

Exotic particles with quantum numbers which can decay into 7.77 are predicted in 
several models [229]. Using a 4.7fb~! data sample at 7, 8 and 13 TeV, the LHCb collabo- 
ration has performed a Dalitz-plot analysis of B? —> n.Kt+m~ decays, where the ne meson 
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is reconstructed in the 7.— pp final state [230]. Evidence was found for a new exotic 
resonance in the ner~ system, later dubbed the X(4100)~ by the PDG. The significance 
of this new resonance exceeds three standard deviations. 


5.7 Structures in the J/yJ/W mass spectrum 


The production of J/)J/w pairs in high-energy pp collisions was observed for the first time 
by the LHCb experiment using a 37.5 pb ' data sample collected in 2010 at \/s = 7 TeV 
231|. The J/J/) mass spectrum was studied in a sample of 116 + 16 signal pairs, and 
no structures was found. The subsequent analysis of 279 pb~! of data collected in 2015 
at ,/s = 13 TeV showed the dominant role of the double-parton scattering (DPS) 
mechanism for J/wWJ/w production over the single-parton scattering mechanism (SPS). This 
in turn includes both a non-resonant SPS contribution and cccé tetraquark production. 
Using the full Run 1 and 2 data-set, the J/J/) mass spectrum was studied in more 
detail [233]. The data, shown in Fig. were found to be inconsistent with the hy- 
pothesis of non-resonant SPS plus DPS in the 6.2 < my jy < 7.4 GeV/c? range, where 
cccé tetraquarks decaying into J/)J/wW pairs are expected. A narrow peaking structure at 
Mappa © 6.9 GeV/c? matching the lineshape of a resonance, and a broader structure near 
to threshold, were found. 
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Figure 5.19: The J/w-pair mass spectrum with the fit result superimposed. The fit accounts for 
the interference between a resonance and non-resonant SPS contribution [233]. 


The global significances of the broader structure close to threshold or the narrow peak 
around 6.9 GeV/c? (provided that the other structure exists), are determined to be larger 
than 5 standard deviations. The structures are consistent with hadron states made up of 
four charm quarks, alternatively they may also result from near-threshold rescattering 
effects, as the XoXo and YX thresholds sit at 6829.4 MeV/c? and 6925.4 MeV/c’, 
respectively. 
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5.8 Light hadron spectroscopy 


n — 7/ mixing has been studied by LHCb in B? > J/n decays, resulting in four 
observed decays modes, using the Run 1 data-set (234). The 7 and 7’ were identified in the 
decay modes n > nnn, n> ntr n?’ and n — p°y. For decays of B? (B°) mesons, the 
n” mesons are formed from initial s3 (dd) quark pairs, hence the measurement of the ratios 
of branching fractions of these decays allows a precise measurement of the n — 7’ mixing 
angle. It also probes the gluonium component in the 7 meson. 


Excited strange mesons have been studied in the ¢K* system from a full amplitude 
fit of Bt J/po K+ decays using the LHCb Run 1 data-set (1424143). Even though no 
peaking structures are observed in the ¢K* mass distributions. correlations in the decay 
angles reveal a rich spectrum of K** resonances. In addition to the angular information 
contained in the K** and ¢ decays, the J/ decay also helps to probe these resonances, 
as the helicity states of the K*t and J/w mesons originating from the Bt decay must be 
equal. Unlike the earlier scattering experiments investigating K* — øK decays, a good 
sensitivity to states with both natural and unnatural spin-parity combinations is achieved. 

The dominant 1* partial wave has a substantial non-resonant component, and at 
least one resonance that has a significance of 7.60. There is also 20 evidence that this 
structure can be better described with two resonances matching expectations for two 
2P, excitations of the kaon. Also prominent is the 2~ partial wave which contains at 
least one resonance at 5.00 significance. This structure is also better described with 
two resonances at 3.00 significance. Their masses and widths are in good agreement 
with the well-established K2(1770) and K2(1820) states, matching the predictions for 
the two 1D 2 kaon excitations. The 17 partial wave exhibits 8.50 evidence for a resonance 
which matches the K*(1680) state, which was well established in other decay modes, and 
matches expectations for the 13D, kaon excitation. This is the first observation of its 
decay to the ¢K final state. The 27 partial wave has a smaller intensity but provides 
5.40 evidence for a broad structure that is consistent with the K3(1980) state, previously 
observed in other decay modes, and matches expectations for the 2? P, state. The K (1830) 
state (31S) candidate), earlier observed in the oK decay mode in K~p scattering, is 
also confirmed at 3.50 significance. Its mass and width is now properly evaluated with 
uncertainties for the first time. 
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6 Measurements not originally planned in LHCb 


While originally designed to study the production and decay of b and c hadrons, LHCb 
has extended its physics programme to also include other areas, such as physics with jets, 
the production of W and Z bosons, searches for new particles in open mode, and nuclear 
collisions. Selected highlights are summarised below. 


6.1 Production of EW bosons W and Z 


LHCb has measured the production of Z and W bosons inclusively and in association 
with jets, reconstructed in mainly muonic final states, using the data collected at ys = 
8 TeV (236). Also decays to ete [237], TETT and ev [239], have been measured, 
however the muon channel is the most efficient due to the excellent performance of the 
muon system (see Sect. |2.6). The Z — ut u~ decay shows a spectacularly clean signal, as 
shown in Fig. a) |240|. The W — uv channel also manifests in a clear signal, shown 
in Fig. [6.1{b) ag The absolute and differential cross sections, their ratios, and charge 
asymmetries have been measured and compared to theoretical predictions. Figure 
(Left) shows the comparison of W and Z cross section measurements to SM predictions, 
showing good agreement. 
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Figure 6.1: (a) Invariant mass distribution of dimuon pairs in the Z-candidate sample. (b) 
LHCb data compared to QCD, electroweak and heavy flavour background, for positive (left) 
and negative (right) muon pr spectra of W candidates |235]. 


6.2 Jets in LHCb 
Measurements of jets at LHCb address several interesting areas of study: 
e Jet properties and heavy-quark jet tagging; 


e The constraining of proton parton density functions (PDFs) and to probe hard 
QCD in a unique kinematic range. Figure |6.3| shows the domain in the (x,Q?) plane 
covered by the LHCb detector, complementing the kinematic ranges of ATLAS and 
CMS; 
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Figure 6.2: (Left) The W versus Z measured production cross sections, showing comparison 
with theoretical predictions. (Right) Contributions to the selected W plus jet sample in the 
discriminating variable ph / pha et (described in the text). 


e Direct searches for the Higgs boson decaying to bb and cē final states; 
e Direct searches for long-lived beyond-the-SM particles decaying into jets. 


Jets are reconstructed in LHCb using a particle flow algorithm clustered using the 
anti-kT algorithm with R = 0.5 |242]. The calibration of jet reconstruction is performed 
in data using Z + yt decays which also contain a jet, where the jet is reconstructed 
back-to-back with respect to the Z. The efficiency for reconstructing and identifying jets 
is around 90% for jets with transverse momentum pr > 20 GeV/c. Furthermore, LHCb 
has developed a method to tag jets and to determine whether they correspond to 
a b or c quark or to a lighter quark. Jets are tagged whenever a secondary vertex (SV) 
is reconstructed close enough to the jet in terms of R =,/(Ad¢? + An?). This provides a 
light-jet mistag rate below 1%, with an efficiency for b (c) jets of ~ 65% (~ 25%). 

Moreover, using the SV and jet properties, two boosted decision trees (BDTs) have 
been developed, one to separate heavy from light jets and one to separate b from c jets. A 
summary of the obtained performance is shown in Fig. where the efficiency of flavour 
identification is plotted as a function of the misidentification of light jets. 


Production of W and Z with jets. W and Z production have also been studied in 
association with jets [236], in W +j, Z +j, W +bb and W + c. Jets are reconstructed 
as described above, whilst Z and W bosons are reconstructed mainly in muonic final 
states. The production of W boson plus jets is discriminated from misidentified QCD 
background processes using a muon isolation variable, which is built as the ratio between 
the pr of the jet containing the muon and the pr of the muon alone. Figure [6.2(b) shows 
the distribution of this variable, with genuine muons from the W boson peaking at 1. 
Figure shows the comparison of the measured cross-sections in LHCb with theoretical 
expectations, showing very good agreement. 

LHCb has measured the WF + bb, W* + cé, production cross sections using a sample 
of pp collisions taken at \/s = 8 TeV with a high-pr isolated lepton from the W decay 
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Figure 6.3: The region in the x — Q? plane probed by LHCb, compared to ATLAS, CMS and 
previous experiments. 


(electron or muon) and two heavy flavour (b or c) tagged jets in the final state. The 
heavy-quark tagging uses the method described above. In this analysis, the W + cce channel 
is studied for the first time. 

In order to extract the different signal components, a simultaneous four-dimensional 
fit is performed on the yt, u~, et and e~ samples. Here the electron channels are used to 
increase statistics. The four variables used in the fit are the dijet mass, a multivariate 
discriminator to separate tt from W +bb and W +cZ events and a multivariate discriminator 
to separate b- and c-jets, used for both accompanying jets [243]. In this fit, the background 
from QCD multi-jets is extrapolated from a control sample in data, while other background 
contributions are fixed to SM theoretical expectations. Only the signal components are 
then unconstrained. The projections of the resulting fit on four input variables for the pt 
sample are illustrated in Fig. The statistical significance of the measured W+ + bb, 
Wt + ce, W~ + bb, W~ + c and tt production cross sections are 7.10, 4.70, 5.60, 2.50 
and 4.90, respectively. The cross sections measured in the LHCb fiducial acceptance agree 
well with the Next-to-Leading-Order (NLO) theory predictions. 
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Figure 6.4: Simulated efficiencies for SV-tagging a b, c-jet as a function of mistag probability for 
a light-quark jet. 


Search for long-lived new particles. The LHCb detector has been designed to 
measure very rare decays of b quarks, with the aim of detecting the presence of new 
beyond-the-SM particles through their couplings in loops, which could change the expected 
SM branching ratios. This implies reaching mass scales higher than those explored in 
open new particle searches, where the particle is produced directly in pp collisions. The 
traditional way to search for new particles is, as for ATLAS and CMS, reconstructing 
their decays through exclusive final states in their invariant masses, or with missing 
energy techniques. Here hermiticity is a mandatory feature of the detector. In all cases, 
assumptions on their coupling, their production mechanism, and on their decay modes 
must be made, and the search is therefore guided by theoretical models. 

For LHCb, the most promising final states are those decays which form a secondary 
vertex, for which the LHCb VELO (see Sect. is extremely efficient, i.e. the new 
particles are long lived. Several negative results have been published which are often 
less sensitive than the General-Purpose Detector (GPD) results. However in some cases, 
LHCb can extend the exclusion region. An example is given in Fig. |6.7| for a Higgs-like 
particle H? decaying into jets forming a separate secondary vertex. The LHCb exclusion 
region is compared with that of ATLAS and CMS, which demonstrates the complementary 
of the LHC experiments. The limits are, in this specific case, competitive, despite a factor 
10 less luminosity. 


Production of tt pairs. Top quark production is an excellent example where the 
forward acceptance of the LHCb detector has several advantages with respect to the 
central region instrumented by ATLAS and CMS. The t quark cross-section can provide 
important constraints on the large-z gluon PDF, where the forward kinematic region 
is particularly sensitive. In addition the forward region provides a greater fraction of 
events with quark-initiated production than in the central region, and enhances the size 
of tt asymmetries visible at LHCb. The challenge for LHCb to measure tt production 
is the small acceptance and the impossibility of a missing energy measurement. Also 
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Figure 6.5: Comparison with theoretical calculations of measured cross sections for W and Z 
production accompanied by jets. The orange bands represent the statistical uncertainty only, 
the yellow bands are the quadratic sum of statistical and systematic uncertainties. 


the fact that the luminosity is limited by the need to reduce multiple interactions for 
measurements in the b sector, disfavours tt statistics. 

Top-quark production is presented here at \/s=13 TeV, which gives an increase in the 
production rate of an order of magnitude with respect to 8 TeV, and which brings these 
new channels into statistical reach. The tt analysis is based on an integrated luminosity of 
2 fb™t, and with eub measured in the final state. Hence the final state is the decay chain 
tt > bWtbW- — et bb, where at least one b-jet is reconstructed. This is a very pure 
final state, as the second lepton suppresses W + bb production and the different flavoured 
leptons suppress Z + bb. The signal purity is illustrated in Fig. [6.8{a), and the observable 
cross section is measured to be g = 126 + 19( (stat)) + 16( (syst)) + 5(lumi) fb, which is 
compatible with SM predictions. 


Z — bb decay. This measurement is an important validation of the LHCb jet re- 
construction and b-tagging performance. Two b-tagged jets are reconstructed, with a 
third balancing jet also reconstructed to help control the QCD background and define 
signal and control regions using a multivariate technique. The background-subtracted 
signal distribution is shown in Fig. [6.8{b) [244]. The signal is observed with a statistical 
significance of 60 and the measured cross section is found to be compatible with SM 
predictions at next-to-leading order. 
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Figure 6.6: Projections of the simultaneous four-dimensional fit for the y+ sample to: 
(a) the dijet mass, (b) the discriminator to separate tt from W + bb and W + cé, and (c) the 
discriminator to separate b and c leading jets (d) sub-leading jets. In light blue is W + bb, in 
green tt, in red W + cé and in black the background. 


6.3 Dark Photons 


The possibility that dark matter particles may interact via unknown forces, almost not felt 
by SM particles, has motivated substantial effort to search for dark-sector forces (see [245 
for a review). A dark-force scenario involves a massive dark photon, A’. In the minimal 
model, the dark photon does not couple directly to charged SM particles, but it can gain a 
weak coupling to the SM electromagnetic current via kinetic mixing. The strength of this 
coupling is suppressed by a factor € with respect to the SM photon. If the kinetic mixing 
arises from processes whose amplitudes involve one or two loops containing high-mass 
particles, perhaps even at the Planck scale, then 1071? < e? < 10~ is expected [245]. 
Constraints have been placed on visible A’ decays by previous beam-dump, fixed- 
target, collider and rare meson decay experiments; the few-loop region is ruled out for dark 
photon masses m(A’) ~ 10 MeV/¢?. Additionally, the region €? < 5 x 1077 is excluded 
for m(A’) < 10.2 GeV/c?, along with about half of the remaining few-loop region below 
the dimuon threshold. Many ideas have been proposed to further explore the [m(A’), e°] 
parameter space, including an inclusive search for A’ + u` wt decays with the LHCb 
experiment. A dark photon produced in proton-proton collisions via y*—A’ mixing inherits 
the production mechanisms of an off-shell photon with m(7*) = m(A’), therefore both 
the production and decay kinematics of the A’ > utu and y* —> pt” processes are 
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Figure 6.7: Comparison of the LHCb exclusion region for the Branching Ratio of H H? > nyny, 
where my is a long-lived particle decaying to jets. Exclusion regions for ATLAS and CMS are 
also shown. 


identical. 

LHCb has performed searches for both prompt-like and long-lived dark photons 

produced in pp collisions at a centre-of-mass energy of 13 TeV, using A’ —> u u` decays and 
a data sample corresponding to an integrated luminosity of 1.6 fb~! collected during 2016. 
The prompt-like A’ search is performed from near the dimuon threshold up to 70 GeV, 
above which the m(u*py~) spectrum is dominated by the Z boson. The prompt-like 
dimuon spectrum is shown in Fig. 
Three main types of background contribute to the prompt-like A’ search: prompt off-shell 
y* — utp, which is irreducible; resonant decays to +, whose mass peak regions 
are excluded in the search (see Fig. (6.9), and various types of misidentification, which 
are highly suppressed by the stringent muon-identification and prompt-like requirements 
applied in the trigger. 

For the long-lived dark photon search, i.e. with displaced dimuon vertices, the stringent 
criteria applied in the trigger make contamination from prompt muon candidates negligible. 
The long-lived A’ search is restricted to the mass range 214 < m(A’) < 350 MeV/c?, where 
the data sample potentially provides sensitivity. In this case the background composition 
is dominated by photon conversions to utu in the VELO, b-hadron decays where two 
muons are produced in the decay chain, and the low-mass tail from K} —> ntr decays 
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Figure 6.8: (Left) The eub invariant mass for all 44 selected tt candidates, illustrates the 
excellent signal purity. (Right) The background-subtracted dijet mass spectrum showing the 


Z— bb signal [244]. 
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Figure 6.9: The prompt-like u" u™ mass spectrum. 


where both pions are misidentified as muons. 

In the dark-photon searches, no evidence for a signal is found, and 90% CL exclusion 
regions are set on the y — A’ kinetic-mixing strength, shown in Fig. The constraints 
placed on prompt-like dark photons are the most stringent to date for the mass range 
10.6 < m(A’) < 70 GeV/c’, and are comparable to the best existing limits for m(A’) < 0.5 
GeV/c?. The search for long-lived dark photons is the first to achieve sensitivity using a 
displaced-vertex signature. These results demonstrate the unique sensitivity of the LHCb 
experiment to dark photons, even using a data sample collected with a trigger that is 
inefficient for low-mass A’ > +p decays. Using knowledge gained from this analysis, 
the software-trigger efficiency for low-mass dark photons has been significantly improved 
for 2017 data taking. 

In Run 3 to come, the planned increase in luminosity and removal of the hardware- 
trigger stage should increase the number of expected A’ — utu decays in the low-mass 
region by O(100— — 1000) compared to the 2016 data sample. 
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Figure 6.10: Results of the dark photon search. Both prompt-like (top) and displaced (centre) 
exclusions are shown. 


6.4 Nuclear Collisions 


Ultra-relativistic heavy-ion collisions allow the study of the so-called Quark-Gluon Plasma 
(QGP) state of matter, a hot and dense medium of deconfined quarks and gluons where 
heavy quarks are crucial probes. Produced via hard interactions at the early stage of the 
nucleus-nucleus collision, before the QGP formation, heavy quarks experience the entire 
evolution of the QGP. A correct interpretation of these probes requires a full understanding 
of Cold Nuclear Matter (CNM) effects, which are present regardless of the formation of 
the deconfined medium. To disentangle the CNM from genuine QGP effects, heavy-flavour 
production in proton-nucleus collisions is studied. 

The LHCb experiment has collected data of proton-lead (pPb) and lead-lead (PbPb) 
collisions. Since the LHCb detector covers only one direction of the full acceptance, there 
are two distinctive beam configurations for the pPb collisions. In the forward (backward) 
configuration, the proton (lead) beam enters the LHCb detector from the interaction 
point. The proton beam and the lead beam have different energies per nucleon in the 
laboratory frame, hence the nucleon-nucleon centre-of-mass frame is boosted in the proton 
direction with a rapidity, y, shift. This results in the LHCb acceptance for the forward 
configuration as 1.5 < y < 4, and for the backward configuration —5 < y < —2.5. 

In addition, LHCb provides the unique capability at the LHC to collect fixed-target 
collisions utilising the System for Measuring the Overlap with Gas (SMOG) system [247]. 
Originally designed for precise luminosity measurements, SMOG provides the injection 
of a noble gas such as argon or helium inside the primary LHC vacuum around the 
VELO detector with pressure O(10~") mbar, allowing measurements of p-gas and ion-gas 
collisions, and operating LHCb as a fixed target experiment. Since 2015, LHCb has 
exploited SMOG in physics runs using special fills not devoted to pp physics, with a 
variety of beam (p or Pb) and target configurations. This allows unique production studies 
which are relevant to cosmic ray and heavy-ion physics. 
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The heavy-ion results on heavy-flavour production in pPb, PbPb and fixed-target 
collisions collected by LHCb bring yet more diversity and complementarity into the 
field. Also in this context, the excellent momentum resolution and particle identification 
provided by LHCb are especially suited for measuring heavy quark production. The LHCb 
collaboration joined the other participants into the LHC heavy-ion collider programme 
with a pPb run at 5 TeV in 2013 and with a PbPb run in 2015. Following these pioneering 
data runs, significantly larger data-sets have been successfully recorded. 


Fixed target collisions. LHCb has reported first measurements of heavy-flavour pro- 
duction with the fixed-target mode [48]. J/ production cross-sections and D° yields have 
been measured in pHe collisions at \/s~n = 86.6 GeV and pAr collisions at \/s~w = 110.4 
GeV, over the rapidity range 2 < y < 4.6. The cross-section measurements are made 
for pHe data only, since the luminosity determination is only available for this sample. 
After correction for acceptances, efficiencies and branching fractions, the cross-sections 
are extrapolated to the full phase space. The D° measurement is used to extract the 
cc cross-section. The J/) and cc measurements are compared in Fig. with other 
experiments at different centre-of-mass energies and with theoretical predictions. 

With pHe data, LHCb also measured the antiproton production cross section [249], 
a very interesting direct determination, helping the interpretation of the antiproton 
cosmic-ray flux detected by space experiments [250]. 
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Figure 6.11: J/w (left) and cé (right) cross-section measurements as a function of the centre-of- 
mass energy, compared with other experimental data (black points). The bands correspond to 
fits based on NLO NRQCD calculations for J/p and NLO pQCD calculations for cZ, respectively. 
More details are given in [248]. 


Collider mode. In collider mode, the LHCb experiment has collected proton-lead 
collision data at syn = 5 TeV in 2013 and at 8.16 TeV in 2016. The 2013 data 
sample corresponds to an integrated luminosity of 1.06 + 0.02nb~! for the forward 
and 0.52 + 0.01nb~' for the backward regions, whilst the 2016 data corresponds to 
13.6+ 0.3 nb~! for the forward and 20.8 + 0.5 nb™t for the backward. These data samples 
are used to measure quarkonium and open charm or beauty production. 

Y(nS)-meson production is studied in the decay to two opposite-sign muons 251). The 
measurements include the differential production cross-sections of T(15), Y(2S) states and 
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nuclear modification factors, performed as a function of transverse momentum and rapidity 
in the nucleon-nucleon centre-of-mass frame of the Y(nS) state. Also the production 
cross-sections for the Y(3S) is measured, integrated over phase space, and the production 
ratios between all three Y(nS) states are determined. 

The three states are well identified in both pPb and Pbp configurations as shown in 
Fig. |6.12| The nuclear modification factors are compared with theoretical predictions, 
and suppressions for bottomium in pPb collisions are observed. The LHCb measurements 
improve the understanding of cold nuclear matter effects down to low pr. 


2600 \Syy=8.16 TeV, pPb } \s.=8.16 TeV, Pbp 
$500 ed 
i=) — 900 E YS) 
Sn S500 Background 
= ~ 
$300f- 24001 —Total 
Sooo 3300 
E "200 
100 Soon 
so 10 11 a" 10 11 
M(u*w) [GeV/c?] M(u+uw) [GeV/c] 


Figure 6.12: Invariant-mass distribution of u" u~ pairs from the (left) pPb and (right) Pbp 
samples after trigger and offline selections. 


7 Future prospects 


Over the years 2011 - 2018, both the LHC machine and the LHCb detector performed 
extremely well, providing great improvements with respect to the B Factory measurements, 
in particular the pioneering CP violation measurements (Sect. 3), the observation of the 
rarest beauty meson decays (Sect.4) and the discovery of pentaquarks (Sect. 5). LHCb 
also observed and reported a number of interesting hints of anomalies related to the 
flavour sector, which has generated much theoretical attention, especially relating to rare 
decays and lepton flavour universality. The precision achieved by the experiment is in line 
with prior expectations, as documented in [252], and which demonstrates the remarkable 
understanding of all aspects of the detector. 

To further pursue these exciting results and fully exploit the flavour physics potential 
of the LHC, the LHCb detector required an upgrade, to increase the rate and efficiency 
of data taking beyond the Long Shutdown 2 (LS2). Consequently the LHCb detector is 
now undergoing a major upgrade that is well underway, and will allow the experiment to 
pursue its superb performance into the future. 

At present, the hardware-based trigger limits the amount of data taken each year to a 
maximum of about 2 fb~!. In addition, most of the detector sub-systems would not cope 
with higher luminosity due to either their outdated readout electronics or radiation-induced 
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damage sustained during Run 1 and Run 2 data taking. The initial ideas regarding the 
upgrade were formulated in 2011 [253], and further solidified in 2012 when the Technical 
Design Report was released |254|. Many of the subdetector components are largely 
unchanged in the upgrade, with the exception of a new pixel vertex detector replacing 
the current VELO, the TT stations being replaced by a new silicon micro-strip upstream 
tracker (UT), and the straw outer chambers replaced by a scintillating fibre detector. 
Details of each subdetector upgrade can be found in refs. [254258]. 

The crucial point of the upgrade project is to build a reliable and robust detector 
capable of operating at higher luminosity without compromising the excellent physics 
performance of the current detector. This, in turn, cannot be achieved by redesigning 
the hardware components alone, but has to be augmented by a new innovative and 
flexible trigger system. A critical part of the upgrade strategy is the design of a so-called 
“trigger-less” front-end electronics system capable of reading out the full detector at 40 
MHz, i.e. at the LHC clock frequency. Completely new and novel chips have been designed 
and tested for the pixel sensors the UT and RICH detectors (257. 

The upgraded detector will operate at an instantaneous luminosity of 2 - 1033 cm?s~ 
which allows collection of around 10 fb~! of data per year as a target, also keeping pace 
with Belle II [259], the other major flavour-physics experiment. Figure shows the 
corresponding time-line for LHCb operations over the next decade. 
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Figure 7.1: A time-line showing the operations of the LHC and the HL-LHC over the next 
decade, including long shutdown (LS) periods, as can be estimated today. The operational 
periods of LHCb and Belle II are shown 


To efficiently run at increased luminosity, the present hardware-based trigger will be 
replaced, and events will be selected by the software-based HLT alone. To cope with the 
much higher event rate (typically five proton-proton interactions per beam crossing), a 
flexible software trigger will be employed and coupled with a re-optimized network capable 
of handling a multi-terabyte data stream. The upgraded trigger will process every event 
(the visible rate at LHCb is estimated to reach 30 MHz) using information from every 
sub-detector to enhance its decision and maximize signal efficiencies, especially for the 
hadronic channels. The precision of particle identification and track-quality information 
will be identical than ” offline” and able to reduce the rate down to 20-100 kHz. The new 
trigger strategy will increase the triggering efficiency for the hadronic channels by a factor 
2 to 4 with respect to Run 1 [260], corresponding to an increase of a factor 10 to 20 for 
the hadronic yields. 

Finally, plans for a further future upgrade (called Upgrade IJ) to use the full potential 
of flavour physics during the HL-LHC operation have now started (261][262]. This upgrade 
would require a complete redesign of the detector able to take data at instantaneous 
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luminosities of 2- 10°* cm?s~!, and collect ~ 50 fb~! of data per year, guaranteeing LHCb 
operation beyond 2030. 

The LHCb Upgrades I and II will significantly improve the reach of key physics 
measurements. By way of example, the precision quoted in Section on today’s 
measurement of the CKM angle y, O(5°), will be improved to 1° with the new Upgrade I 
hadronic trigger and luminosity increase. This further improves to 0.35° with the large 
statistics accumulated with LHCb Upgrade II. As discussed in Section [4.1] currently there 
is not enough sensitivity to measure the rare decay B(B° > p*p~). With LHCb Upgrade 
I, a first observation should be possible, but it will require Upgrade II to reach a ~ 10% 
precision on its branching ratio. This sensitivity will allow a meaningful measurement of 
the ratio between B° and B? into the wtp final state, and will constitute a clean and 
powerful test of extensions beyond the SM. 


8 Summary and Conclusions 


In ten years of operation at the LHC, the LHCb experiment has delivered a remarkably 
rich programme of physics measurements. In this paper the 25-year evolution of the 
experiment since its inception has been described, and its successes and achievements 
have been summarised. The diversity of the physics output has truly shown LHCb to be 
a “general-purpose detector in the forward region”. 

Over the last ten years, LHCb has measured the CKM quark mixing matrix elements 
and CP violation parameters to world-leading precision in the b- and c-quark systems. 
The experiment has measured very rare decays of b and c mesons and baryons, some with 
branching ratios down to order 107°, testing Standard Model predictions to unprecedented 
levels. Hints of new physics in rare-decay angular distributions and through tests of 
lepton universality in electron-muon decay modes have generated considerable theoretical 
interest. The global knowledge of b and c quark states has improved significantly, through 
discoveries of many new resonances already anticipated in the quark model, and also by the 
observation of new exotic tetraquark and pentaquark states. In addition, many interesting 
measurements have been made that were not anticipated in the original LHCb proposal, 
such as electroweak physics, jet measurements, new long-lived particle searches and heavy- 
ion physics. An incredibly rich harvest of fundamental results has been produced, many 
of these will remain in textbooks for years to come. 

LHCb has recently been upgraded and will start data-taking early in 2022 at a factor 5 
higher luminosity, incorporating new subdetectors and a software-based trigger. Statistics 
in hadronic modes will be improved by at a factor 10-20, allowing much more precise 
measurements, especially of very rare b- and c-hadron decays. In addition, the future 
planned Upgrade II at the (HL)-LHC early in 2030 will ensure that LHCb maintains its 
lead in flavour physics for at least the next two decades. 
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